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INTRODUCTION 


The 7- by 10-foot high speed tunnel is equipped with a digital data 
acquisition, display, and control system operated by a dedicated on-site 
computer. The computer software consists of three main parts: the real 

time batch monitor operating system (RBM) , the operating acquisition program 
(OAP) , and the real time applications task (RTAT) . Real time data reduction 
is accomplished by RTAT which runs as a task under OAP. 

The 7- by 10- foot high speed tunnel performs a wide range of tests 
employing a variety of model installation methods. To support the reduction 
of static data from this facility, a generalized wind tunnel data reduction 
program had been developed for use on the Langley central computer complex. 
This program, which was developed over a number of years by a data reduction 
support contractor, contains the best algorithms available for wind tunnel 
data reduction. Since RTAT was subject to the same requirements for 
generality and flexibility as the central program, it was decided to derive 
RTAT from the central program instead of writing a complete new program for 
the on-site computer. Thus, RTAT is a generalized wind tunnel data reduction 
program which contains the best algorithms available. The user invokes 
those RTAT features necessary for a particular test through the input 
specifications. To fit into the available memory on the on-site computer, 
RTAT uses a more sophisticated data management scheme than the central 
program and therefore the input specifications are slightly different. 

Tliis report describes the capabilities of the current release of RTAT. 

It is designed as a user's guide, an operator's manual, and a reference 
manual. It is not a programmer's guide. Thus, it contains detailed 
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descriptions of the :nput specifications, instructions for the console 
operator, and full descriptions of the algorithms but does not document 
t’ne code itself. A detailed description of the complete hardware and 
software environment in which RTAT executes is beyond the scope of this report. 

The purpose of this report is to assist personnel involved in the conduct 
of tests in the 7* by 10-foot high speed tunnel. 

HARDWARE 

The 7- by 10-foot high speed tunnel digital data acquisition, display, 
and control system presently consists of a Honeywell (Xerox) Sigma 3 computer 
with 49152 words of memory in the central processing unit, an external input/ 
output processor, 4.5 megeibytes of rapid access disk storage, two nine- 
track tape drives, a card reader, a line printer, a Tektronix 4014 graphics 
terminal and hardcopy unit, a data acquisiton unit, and a data link to the 
central computer complex. 

The data acquisition unit presently supports 50 analog channels, 20 digital 
channels, and 8 tachometer channels associated with the tunnel and 50 analog 
channels and 10 digital channels associated with a static calibration area. 

A number of similar systems exist at other Langley research facilities. 

SOFTWARE 

The operating system used on the Sigma 3 is the Honeywell (Xerox) 
real time batch monitor (RBM) system, A detailed description of RBM 
is beyond the scope of this report. 

The operating acquisition program (OAP) was obtained as part of a 
combined hardware/software procurement for a number of similar data acquisition 
systems. The OAP runs as a foreground resident program under RBM. A 



brief description of GAP operation and GAP input specifications is 
provided in APPENDIX A. RBM and GAP provide a fixed environment in which 
RTAT must execute. 

RTAT was developed specifically for the 7- by 10-foot high speed 
t’onnel. RTAT consists of a frame task and a cyclic task which are written 
in FORTRAN and are assembled as part of OAP. 

The RTAT frame task performs a fast calculation of Mach number and 
dynamic pressure. These vlaues are displayed at the tunnel drive control 
operator’s console and updated approximately t^'n times per second. 

The RTAT cyclic task performs all of the data reduction computations. 

The results are displayed at the data acquisition operator's console and 
at the model attitude control operator's console and are updated approximately 
once per second. 

The data base concept employed by RTAT is described in APPENDIX B which 
includes a glossary of standard vari=^ble names used by RTAT. 

The general input card specifications for RTAT are described in APPENDIX C. 
The extra equation input card specifications are described separately in 
APPENDIX D. The input specif ications for interactive calibration workup 
sessions are described in APPENDIX E. 

Some sample RTAT input setup decks are presented in APPENDIX F. 

A description of the computational algorithms is presented in 
APPENDIX G. 

The RTAT operating procedures are discussed in APPENDIX H. 

Some hints on performing a manual check of the computations executed 
by RTAT are given in APPENDIX I. 
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RTAT IDLE LOOP 


The idle loop is the path followed by RTAT to update the real time 
displays between data points. Raw data are obtained from the latest fraune 
sampled by the OAP. These data are fully processed exactly as it would be 
during normal execution with the following exceptions: 

(a) all I/O operations to the card reader, line printer, graphics 
device, and tape drive are suppressed, 

(b) wind-off zero and tare data are not saved as they would be for a 
data point, and 

(c) unless specifically requested in the input specifications, 
second order balance interactions are omitted. 

RTAT NORMAL EXECUTION LOOP 

The normal execution loop is the path followed by RTAT in response to a 
discrete event. Raw data are obtained from the latest OAP averaged record 
buffer. The data are fully processed and appropriate I/O operations 
performed. The following sections will describe the normal execution loop, 
more or less following the functional layout in Figure 1. 

RTAT SPECIAL FUNCTION EXECUTION 

In addition to the routine execution performed within the idle loop 
and normal execution loop, RTAT may bypass those loops and execute a specific 
function in response to a user request. Some examples of specific functions 
are weight tare computations, calibration data acquisition, calibration data 
computations, and system calibrate computations. 
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ktat planned future expansions 


RTAT will support powtiiCd i^ndel testing using high pressure air to supply 
exit nozzles on the model. However, that portion of RTAT intended to support 
the complexities of static powered calibrations and the interactive workup of 
static powered calibration data has not been coded. RTAT will add the code 
to support these functions. 

The hardware capabilities of the 7- by 10-foot high speed tunnel data 
acquisition system are being expanded or are proposed to be enhanced. The 
possible future expansions to RTAT to support these capabilities are: 

(a) A data link between the Sigma 3 and the Langley Central Computer 
Complex is presently being installed which will permit the Sigma 3 to submit 
priority batch data reduction jobs to the central computer complex. RTAT 

will add to the special function execution capability to support this data link. 

(b) The model attitude control system is being converted from manual 
operation to automated operation. RTAT will add to the normal execution 
loop capability to control the automatic operation. 

(c) The mechanically multiplexed pressure transducers may be replaced 
with electronically multiplexed pressure transducers. When that occurs, 

OAP and RTAT will be updated appropriately. 

(d) The entire Sigma 3 data system may be replaced with a current 
state of the art data system. Since RT;^T is entirely written in FORTRAN, 
employs completely general algorithms, and has a relatively straightforward 
interface into the system dependent features of the Sigma 3 system, RTAT 
will be converted to the new data system. After conversion, certain aspects 
of RTAT may be revised to take advantage of additional capabilities in the 
new data system. 
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Figure 1.- RTAT Functional Block Diagram, 
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Figure 1,- Concluded 
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AN OVERVIEW OF OAP 

OAP operates the data acquisition unit hardware and directs the acquisition 
of the raw data. A frame of raw data consists of a single set of channel 
readings. Both OAP and RTAT use the frame data for idle (frame) loop 
computations. Upon the command to acquire a data test point, OAP will assemble 
frames of data and average them appropriately to form a data test point. The 
test point data are used by OAP and RTAT for normal execution (cyclic) 
loop computations which include the generation of output data. 

OAP runs as a resident foreground program and consists of a number of 
tasks attached to different interrupt levels- OAP triggers the frame 
interrupt and the cyclic interrupt to which RTAT is attached. A task will 
begin executing each time the associated interrupt is triggered while the 
interrupt is armed and engibled. jAn interrupt is disabled while the associated 
task is executing and is enabled again when the task is completed, 

OAP includes a permanently resident blank common area which is used for 
communication between tasks. 


OAP Operation 

Load an OAP setup deck into the card reader and start it up. When 
OAP begins execution, a message will be output on the teletype requesting 
that a System Control Panel ENTER be input. At this time, all control 
panels should be placed in the configuration needed for the test. For 
OAP to start, the System Control Panel must be setup with SELECT OPTIONS and 
REQUEST on. An attempt to start OAP without these selections will result in 
an error message output on the typewriter. At this pc.'nt a System Control 
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Panel ENTER should be input. GAP will then read and process the GAP setup 
cards and any errors wilJ be output on the line printer. GAP will read all 
of the control panels and perform a consistency check between the setup cards 
and the hardware configuration. If any discrepancies are detected, the 
errors will be output on the typewriter. Before proceeding further, it is 
necessary to clear all of these errors. This is accomplishea by correcting 
the setup via cards or typewriter and/or changing the panel settings. 
Depressing a panel ENTER button will cause GAP to perform a new consistency 
check . 


Analog System Calibration 

There are two system calibration functions that can be performed 
by the operator: 

1. Normal Calibrate 

2. Hot Zero or Offset Calibration 

A normal calibration is obtained by selecting CALIBRATE on the 
System Control Panel, setting the DATA IDENT thumbwheels to 00, and then 
depressing the System Control Panel ENTER button. The GAP will then 
calibrate all active analog channels and store the results on RAD file 
CALDATA as well as reporting the results on the line printer. 

A hot zero or offset calibration is obtained in the same way as a normal 
calibration but with a DATA IDENT thumbwheel setting of 03. This calibration 
results in GAP storing an offset for each active channel which is subsequently 
applied to every frame of data. Obviously this is only valid if all analog 
channels contain differential transducers. Since the 7- by 10-foot high speed 
tunnel normally uses a number of absolute analog transducers, this type of 
calibration should not be used. 
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If during the calibration^ the operator should desire to cancel the 
calibration, this may be accomplished by depressing the RESET and ENTER 
buttons on the System Control Panel, This will cause OAP to abort the 
cc'’ ibration. 


Test Point Recording 

To cause OAP to take a test po:nt, the operator must select DAU DATA, 
set the DATA IDENT thumbwheel value to any value between 00 and 96, depress 
the ENTER button on the System Control Panel, and tneii depress the DAS 
READY/CYCLING button for each test point desired. 

Test Point Intervention 

If during the recording of a test point, the operator should desire to 
terminate the test point, this may be accomplished by depressing the 
DAS READY/CYCLING button again. 

It is aJso possible to correct the OAP setup during a test point by the 
following steps: 

(a) Depress the HOLD/CONTINUE button on the System Control Panel to 
place the system in HOLD. 

(b) Correct the setup 

(c) Depress the RESET and ENTER buttons on the System Control Panel 

(d) Depress the HOLD/CONTINUE button on the System Control Panel to 
place the system in CONTINUE. 

Selecting LIST on the System Control Panel will cause OAP output to go 
to either the line printer or the typewriter as specified in the OAP setup* 
Selecting TAPE on the System Control Panel will cause OAP to output data for 
a test point to the tape. Setting a channel number in a display thumbwheel 
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and pressing the associated ENTER button will cause GAP to display the 
value associated with that channel. 

GAP Input Card Specifications 

The GAP setup cards follow a very restricted free format. The setup 
cards consist of a card type identifier in column 1 and one or more fields. 
The card type identifier and the fields must each be separated from the nart 
by one or more blanks. More than one field is legal, but fields must be 
separated by one or more blanks. Column 72 must contain either a blank or a 
right parenthesis. Columns 73 through 80 may contain any information the 
user desires. Each field consists of two pcirts: the FIELD IDENTIFIER 

and the FIELD ARGUMENT. Multiple arguments are permitted within a single 
field, but the field must then be enclosed in parenthesis and the arguments 
separated by commas. Note that blemks may not be embedded in a field and 
that an argument may not end with a decimal point. 

The control card type identifiers and the card description are given 

Description 
INTERRUPT CGNTROL 
CHANNEL PARAMETER 
DISPLAY SELECTION 
CONTROL PARAMETER 
DIGITAL CONSTANT SPECIFICATION 
APPLICATION PROGRAM SPECIFICATION 
INPUT DATA DELIMITER 


below: 

Type 

1 
1 

2 
3 

5 

6 
9 


n 
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ThtL only setup cards required by OAP are a type 1 card and a type 9 card. 

In the absence of any other card types or of any field, OAP will use the 
default specified in the descripcien of the card. 

A description of all card types is given below which employs the 
follovjing conventions: 

(1) a - represents an alphabetic character in a field entry 

(2) m or n - represents a numeric digit in a field entry. 

The plus sign, minus sign, and decimal point are acceptable in a numeric field. 

(3) X “ represents any character in a field entry. 

(4) The field arguments are represented at the maximum length 
appropriate to each parameter by the number of characters shown for each 
field argument. An actual argument may use fewer characters. 

Setup Card Type I, Interrupt Control 
This type can only be input by a card. The interrupt card contains 
two f-i-elds, the first gives the interrupt level ncime and the second gives 
the interrupt level to which that name is to be assigned. One and only 
one interrupt level is allowed per card. 

The OAP interrupt defaults listed on the following page are usually 
sufficient. 
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OAP Task Interrupt Assiqrimonts 


NAME 


LEVEL 


DATCHN 0 
rXVICH I 


CTL PNl.S 

ecu 

DAUCTRl.T 

LAMPCTKL 

inSPPNLT 

DISPPCVT 

MSTSDAIT 

MSTSBKUP 

MSTSETUP 

EKKTRAPT 

CAPUNSTO 

OAP’nVDTC 

SYi^tTKLT 

OTSTCTRL 

MSTSDCOI 

DICCONSl 

DIc;CONSJ 

TSl FRAME 

TSIU’RAMK 

TSICYCLT 

TS-VYCLT 


OAP 07 
RBM HO 


110 
111 
IIJ * 

113 

114 

115 
11b 
117 
11 3 
110 
llA * 
IIB 
lie 
iin 
HE 
HE 
120 

1 21 
122 * 

123 * 

124 * 

125 
12b * 
127 
12H * 
129 * 
12A * 
12B * 
12C * 
12D ♦ 
12E 
12F 


Task Assiqnmonts Valid 
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Setup Card Type 1, Channel Parameters 


CHNO , nnn 

This entry is mandatory as the first field on all type 1 cards. The 
channel numbers are assigned to the different channel types as follows: 
Channel Numbers Type 


1 - 100 
151 - 190 
211 - 222 


Analog Channels 
Digital Channels 
Tachometer Channels 


Only one channel number may be specified per card. 


NAME , xxxxxx 

User may assign names to specific channels for use in printer and/or 
typewriter displays and in referring to the channel externally; however ^ only 
a total of 30 names are allowed. 

DEFAULT OPTION: None 


OFSTr nnnnnn 

This permits the user to override the calibration calculated offset. 

The value argument is given in millivolts, the allowable values ranging from 
a maximvun of 65.535 to a minimum of -65.536. A decimal point is required. 
DEFAULT OPTION: Calibration offset is used. 
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UPLM, nnnnnn 

Upper limit for alarm condition check. For analog channels, it is 
stated in millivolts, the legal values ranging from a maximum of 65.535 to a 
minimum of “65.536, For the digital channels and the tachometer channels, 
the value is stated in counts, the legal range being from zero to 99999, 

Only a total of 30 limits are allowed. A decimal point is required- 
DEFAULT OPTION: No check is made. 


LWLM, nnnnnn 

Lower limit for alarm condition check. For analog channels, it is 
stated in millivolts, the legal values ranging from a maximum of 65.535 to a 
minimum of -65.536, For the digital channels and the tachometer channels, 
the value is stated in counts, the legal range being from zero to 99999. 

Only a total of 30 limits are allowed. A decimal point is required. 

DEFAULT OPTION: No check is made. 

DLTE, n 

To delete limit checking set n = 1. 

DEFAULT OPTION; Limits are checked. 

SORT, nnnnn 

Scan rate for continuous scanning patterns. It is stated in SAMPLES/ 
SECOND, the range being from 1 to 20,000. Although this is not presently 
implemented, this field identifier is accepted by OAP and its presence is not 
considered an error. 


DEFAUI T OPTION : None . 
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RNGE, XX 

Specifies whether the ranging should be a set range or automatic. 

Valid parameters for this are A, 8, 16, 32 and 64. If A is specified, 
automatic ranging is used. 

DEFAULT OPTION; omatic rang, 'ig (A) is used. 

PVID, X 

Pressure valve indicator ajisociat js specific hardware with the channel. 

An entry of 0 (zero) causes the relaU i channel to no longer be treated 
as a pressure valve channel. Up to 6 channels may be designated as 
pressure valve channels. 

DEFAULT OPTION: 0 (zero), non pressure valve channel. 

FSET, n 

This is used to specify the det ired setting of the filter associated 
with that analog channel. The valid parameters and the filters which they 
specify are; 

n_ filter 

1 TYPE 1 

2 TYPE 2 

S SPARE 

W WIDE BAND 

It must be noted that the system wi"^ : not proceed past the setup phase until 
all the hcurdware settings of the fil^ 3rs are in agreement with the settings 
specified by the input pareuneters. 

DEFAULT OPTION: W (WIDE BAND) 

‘6 
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LCE (nnnnnnn, nnnnnnn) 

Linear conversion coefficients. Two numeric arguments are required 
for this entry. Only one pair of linear coefficients may be specified per 
channel. Only 30 linear conversion coefficients are allowed, 

DEFAULT OPTION ; None 

TSET, n 

This is used to specify the time period in seconds setting for the 
tachometer channels. The valid parameters are: 

E period 

1 0.1 
2 1.0 
3 10.0 

DEFAULT OPTION; 1 

Setup Card Type 2, Display Selections 

LLIM, nnn 

Listing limit for Type 2 calibration listing. Valid parameters range 
from 0.01 to 0,99, 

DEFAUI.T OPTION; Type 1 calibration listing is used. 

PLUS , nnn 

Percentage of full scale to be used in the calibration routine for all 
four ranges. It is specified as a decimal number, the legitimate range being 
from 0.1 to 0.99. 

DEFAULT OPTION: 0.90 
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SCOD, a 

System configuration, channel data and calibration output destination. 
Either P (Printer) or T (Typewriter) may be designated. 

DEFAULT OPTION: P (Printer) 

Setup Card Type 3, Control Parameters 

RTRY, nn 

Number of retries to be attempted if parity error occurs in MSTS data 
transmission. From 1 through 99 may be specified. MSTS data transmission 
is not implemented. 

DEFAULT OPTION: 7 

AVG, nn 

Length of averaging period, stated in seconds, if scheduled averaging is 
used. Valid parameters range from 1 to 65. 

DEFAULT OPTION: None. Required input for scheduled averaging. 

SCAN, nnn 

Scan rate for discrete scanning pattern. Stated in frames/second, 
from 1 to 500 being acceptable. Entries evenly divisible into 2000 are 
recommended. 

DEFAULT OPTION: 10 frames/second 

PSSR, nnnnn 

Pressure valve stepping rate. Stated in seconds, this is the time 
interval to be allowed for valve stepping and settling in programmed stepping 
rate. Valid entries are numbers from 0.2 through 100.0. 

DEFAULT OPTION: None. Required input for programmed stepping rate. 
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PVCY, nnn 

This indicates the number of times to cycle through the scani- 
valves, recording one sample-per-port; nnn may be any integer from 0 
through 533. An entry of zero (0) causes this to be treated as a null 
field and the operation becomes normal. 

DEFAULT OPTION; 0 (zero) 

TPSC, X 

This entry causes each frame of data to be considered a test point. 

Any enrry is legal, but an entry of zero (0) causes this to be treated as a 
null field and operation returns to normal. 

DEFAULT OPTION: 0 (zero) 

An (mm, nnnnnnn, ) 

This entry builds tables of angles for use by application programs; 
n in the name field identifies the table and must be from 1 to 6. The 
field arguments are in pairs, the first entry giving the position in the 
table, the second entry giving the value of the entry. Thus, (mm, nnnnnnn), 
where mm is any value from 1 through 32, and nnnnnnn is any floating point 
number up to 7 digits. This card allows storage of up to 6 tables with 
32 floating point values each, i.cr each test side. 
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Setup Card Type 5, Digital Constant Specification 

Type 5 cards assign named variables to specific thumbwheel switches on 
the digital constants panel. A maximum of 20 type 5 cards may be included 
in the setup deck 

ILOC, xxxxxx 

xxxxxx is the name given to the constant; i.e., the internal location 
of the constant. 

CON (nn,nn) 

The two field arguments specify the switch numbers of the first and 
last switches of a consecutive set of switches on the DIGITAL CONSTANTS PANEL 
whose settings will determine the value assigned to the constant. The 
value must be no less than 9 and no greater than 48, with the first value 
being less than or equal to the second value. 

Setup Card Type 6, Application Program Specification 

TIME, nnnn 

nnn is any number from 0.01 through 100.0 and gives the cycle time 
interval in seconds for triggering of the cyclic application interrupt. 

DATE, x 

Deletes the cyclic program. 

FRM, n 

For n = 1, causes OAP to trigger trame tasks. 
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Setup Card Type 9, Input Data Delimiter 

END 

This card notifies GAP that setup is complete. 

Use of Typewriter for OA? Setup Input 

With the exception of setup card Type 1 , all other legal setup entries 
may also be made by means of the console typewriter. Typewriter input is 
limited to 72 characters per line. In case of duplication of entries, the 
last entry is the one recognized by GAP. 

Additional Considerations 

The setup deck must be followed by a card, preferably one containing 
!egD in columns 1 through 4, although this is not mandatory, even a blank 
card being acceptable to the card read routine. 

Although GAP does check for certain logical inconsistencies, it accepts 
several versions of the same card, storing one set of values in over the 
other; the last value input being the one used. This allows typewriter 
input to override card input at any time. 

One of the places where GAP does refuse to accept inconsistency is 
the FSET parameter. Until the filter settings on the channels correspond 
to the settings specified, the program will not proceed. Either the filters 
must be reset or the FSET cards must be changed (by typewriter input) . 
Similar considerations apply to the TSET parameter. 

As the setup cards are read, a listing will be printed of those cards 
containing errors. Cards which are correct will NOT be printed, but faulty 
ones will be shown with the errors flagged with dollar signs. The nature of 
the mistakes will be given on subsequent lines. Fault arguments will be 
ignored, but may be corrected through the typewriter input. 
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When specifying arguments given in millivolts or counts, certain routines 
truncate and cause a loss of accuracy when numbers are converted to floating 
point. Therefore, numbers of small magnitude (e.g., 0.002) cannot be 
relied upon to act as triggers or limiting values. 
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RTAT DATA BASE CONCEPT AND GLOSSARY OF NAMES 
Since RTAT was derived from the central computer complex data reduction 
program, it is necessary to understand certain fundamental concepts which 
are employed in the central program - 

The data reduction process is essentially a sequential computational 
process which may involve the use of a number of different programs. 

The standard interface file (SIF) concept was originally developed as a 
standard intermediate tape format to provide communication between diverse 
programs. On a SIF all records have the format: 

XKEY,N, (DATA(I) , I = 1,N) 

where XKEY indicates the type of record. The initial record must be a name 
record; subsequent records may be any type including name records. For any 
non-name record type and for 1 £ K _< N, the contents of DATA(K) will be the 
value asso iated with the ncime contained in DATA(K) in the most recently 
preceding name record. The SIF concept defined a sequential name orir^ited 
file structure providing variable length records having arbitrary internal 
organization. The SIF concept essentially defined an external data base 
which came to be viewed as a natural data base underlying the data reduction 
process . 

As new programs such as the central data reduction program were written, 
they adopted the SIF concept as the internal data base driving the code. 

Each module would obtain its inputs from the data base and store its 


answers in the data base. 
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RTAT uses the same basic SIF concept for both an external and internal 
data base. However, the limited memory available on the Sigma 3 dictated 
several modifications to the internal data base SIF concept described 
above. KTAT incorporates certain keywords as non-SIF constants in all 
record types in the internal data base. Instead of a static data base, 

RTAT employs a dynamic implementation of the data base in which the input 
specifications and the data base share the same area in memory. 

RTAT maintains a copy of the input specifications on a disk file, 
maintains a names record on a disk file, and maintains a non-naune or values 
record on a disk file. The first execution of RTAT for a new setup generates 
the names record and the access map. Each module searches the names record 
fot its inputs, maps their locations, and reserves locations in the data 
base for its answers. RTAT maintains the access map in a permanently 
resident COMMON area. On subsequent executions, modules use the map to 
access the values record. The input specifications, names record, and 
values record share a permanently resident COMMON data base area. At the 
start of each execution, RTAT copies the input specifications into the data 
base area and then appends the raw data obtained from OAP. RTAT then 
t ^»cutes the computational algorithms. Since the input specifications are 
ordered ir the sequence in which they are used, RTAT knows when it ^las 
finished using a block of input specifications. RTAT then collapses or 
shifts that block of input specifications out of the data base area thereby 
making space available for the results which are being appended to the data 
base. When RTAT completes executing the algorithms, all input specifications 
have been collapsed out of the data base area and only raw data and results 
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remain in the data base which can then be written on a SIF answer tape. When 
the next execution of RTAT begins, RTAT refreshes the data base by copying 
the input specifications into the data base area and then entire process is 
repeated. 

The following glossary defines all of the standard names used by RTAT. 


Standard Name 


Description 


ms 

ACOS 

AF 

AF2 

AFB 


AFBA 

AFC 

AFCH 


AFMA 

AFTA 

AFTO 

AFWA 


Extra equation type identifier 
Extra equation type identifier 
Uncorrected delta axial force component 
(Main Balance) 

Uncorrected delta axial force component 
(Second Balance) 

Axial force in pounds due co base pressures 
which is to be applied as a correction to 
balance loads 

Balance axis axl<il ^orce component 
Correct delta axial force component 
Axial force in pounds due to chamber pressures 
which is to be applied as a correction to 
balance loads 

Model axis axial force component 
Axial force weight tare component 
Correct total axial force con^nent 
Vlind axis axial force 
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Standard Name 

ALPG 

ALP I 

ALPS 

ALPU 


ALPW 

ALPZ 


ASIN 

ASQn 

ATAN 

ATN2 

B 

BALA 

BETA 

BETn 

BLK 

BTCH 

CA 


Description 

Gravity axis pitch angle in degrees 
Initial sting pitch deflection angle in degrees 
Sting pitch deflection angle in degrees due 
to load 

Input D-type constant giving the upflow angle 
in degrees as a pitch rotation from the 
wind axis to the gravity axis 
Angle of attack in degrees 
Pitch attitude of the balance with respect 
to gravity when the wind-off zero is 
recorded, in degrees 
Extra equation type identifier 
Cross-sectional area of throat of Flowmeter n 
Extra equation type identifier 
Extra equation type identifier 
Input L -type constant giving model reference 
span in inches 

Input keyword. A balance name and calibration 
date immediately follow this keyword. 

Angle of sideslip in degrees 
Diameter ratio of throat to inlet of 
Flowmeter n 

Input keyword controlling blockage and jet 
boundary correction 
Batch number 

Model axis axial force coefficient 








standard Name 

CAI 

CAB 

CAC 

CAS 

CBAR 

CD 

CDPR 

CDW 

CL 

CLC 

CLSQ 

CLW 

CM 

CMS 

CMW 

CN 

CNI 

CONV 

COS 

CP 


APPENDIX B 

Description 

Extra equation type identifier 
Coefficient of axial force due to chamber 
pressures which contributed to AFB 
Coefficient of axxal force due to chamber 
pressures which contributed to AFCH 
Extra equation type identifier 
Input D-type constant giving model reference 
chord in inches 

Stability axis drag coefficient 
Blockage drag factor 
Wind axis drag coefficient 
Stability axis lift coefficient 
Blockage lift factor 

Stability axis lift coefficient squared 
Wind axis lift coefficient 
Model axis pitching moment coefficient 
Stability axis pitching moment coefficient 
Wind axis pitching moment coefficient 
Model axis normal force coefficient 
Extra equation type identifier 
Input keyv^ord. If nonzero, balance 

interaction iterations will be printed out. 
Extra equation type identifier 
Extra equation type identifier 
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Standard Name 

CPBn 

CPCn 

CRM 

CRMS 

CRMW 

CSS 

CY 

CYM 

CYMS 

CYMW 

CYS 

CYW 

Din 

D2n 

DCn 

DELA 

DELM 

DRAG 

ELSE 

ENDD 


ENDF 


Description 

Pressure coefficient for base pressure n 
Pressure coefficient for chamber pressure n 
Model axis rolling moment coefficient 
Stability axis rolling moment coefficient 
Wind axis rolling moment coefficient 
Extra equation type identifier 
Model axis side force coefficient 
Model axis yawing moment coefficient 
Stability axis yawing moment coefficient 
Wind axis yawing moment coefficient 
Stability axis side force coefficient 
Wind axis side force coefficient 
Inlet diameter of Flowmeter n 
Throat diameter of Flowmeter n 
Discharge coefficient of Flowmeter n 
Jet boundary angle of attack correction 
Jet boundary pitching moment correction 
Stability axis drag 
Extra equation type identifier 
Input optional keyword following real time 
display specifications. If present, all 
constants preceding ENDD will be collapsed 
out of the data base. 

Input optional Keyword following the balance 
interaction deck. If present, all constants 
preceding ENDF will be collapsed out of the 
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standard Name 


ENDF (concluded) 


ENDG 


ENDH 


ENDM 


ENDP 


ENDQ 


ENDS 


APPENDIX B 

Description 

data base, ENDF is required after each 
interaction deck if NBAL > 1. 

Input required Keyword following the summary 
print specifications- All constants 
preceding ENDG will be collapsed out of 
the data base. 

Input required Keyword following reorder 
hookup specifications. All constants 
preceding ENDH will be collapsed out of 
the data base. 

Input required keyword following special 
point-by-point plot specifications. All 
constants preceding ENDM will be collapsed 
out of the data base. 

Input required Keyword following the point- 
by-point print specifications. All 
constants preceding ENDP will be collapsed 
out of the data base. 

Input optional Keyword following the engineering 
unit specifications. If present, all 
consLants preceding ENDQ will be collapsed 
out of the data base. 

Input required Keyword following summary 
plot specifications. All constants 
preceding ENDS will be collapsed out of the 
data base. 
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Standard Name 


ENDT 


ENDU 


ENDX 


EXP 

Fn 

FDPn 


FPn 


FP2n 

FTn 

GPn 


HI 

ID 


Description 

Input required Keyword following normal 
point-by-point plot specifications. All 
constants preceding ENDT will be 
collapsed out of the data base. 

Input required Keyword following the engineering 
unit extra equation specifications. All 
constants preceding ENDU will be collapsed 
out o^ th^- data base. 

Input required Keyword following the force 
extra equation specifications. All 
constants preceding ENDX will be collapsed 
out of the data base. 

Extra equation type identifier 

Velocity of approach factor of Flowmeter n 

Flowmeter differential pressure in psi of 
Flowmeter n 

Flowmeter inlet static pressure in psi of 
Flowmeter n 

Throat static pressure of Flowmeter n 

Flowmeter temperature in of Flowmeter n 

Input Keyword associated with summary print 
specifications 

indicated tunnel total pressure in PSF 

Data identification code number 
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Standard Name 

IFEQ 

IFGE 

IFGT 

IFLE 

it^r 

INTR 


lORD 


J2 

J3 

K 

KB 

KBI 

KDFL 


Description 

Extra equation type identifier 
Extra equation type identifier 
Extra equation type identifier 
Extra equation type identifier 
Extra equation type identifier 
Extra equation type identifier 
Input Keyword signifying that an I-card 

format balance interaction deck immediately 
follows 

Input optional Keyword giving the order of 
interactions to be considered for the 
balance. If omitted^ second order is 
assumed for data points and first order 
for the idle loop. 

Input D-type constant specifying angle of 
attack jet boundary correction factor 
Input D-type constant specifying pitching 
moment jet boundary correction factor 
Total blockage correction factor 
Body blockage correction factor 
Input D-type constant specifying body 
blockage correction factor 
Input D-type constant specifying the 
deflections are to be calculated on the 
basis of total loads or delta loads. If 
omitted, TOTAL is assumed. 
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Standard Name 

KW 

KWI 

KWK 

L/D 

LAMn 

LIFT 

MACH 

MACH 

MPR 

MPLT 


NBAL 


NBAS 


NBM 


NCBR 


Description 

Wing blockage correction factor 
Input D-type constant specifying wing 
blockage correction factor 
Wake blockage correction factor 

Stability axis lift to drag ratio 
Normalized pressure of Flowmeter n 
Stability axis lift 
Free stream test section Mach number 
Extra equation type identifier 
Uncorrected Mach number 

Input Keyword giving the number of special 
point-by-point plot specifications which 
must immediately follow. The maximum is 4. 
Input Keyword giving the number of balances 
to be computed. The current maximum is 2. 
Input Keyword giving the munber of base 
pressure coefficient specifications which 
must immediately follow 
Input Keyword giving the number of rotation 
specifications from the balance axis to the 
model axis which must immediately follow. 
The maximum is 12. 

Input Keyword giving the number of chamber 
pressure coefficient specifications which 
must immediately follow. 
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Standard 

NCP 

NDSP 

NEU 

NEXF 

NEXU 

NF 

NF2 

NFB 

NFBA 

NFC 


N ame Description 

Input Keyword giving the number of pressure 
coefficient array specifications which 
must immediately follow 

Input Keyword giving the number of real time 
display specifications which must 
immediately follow. The maximum is 100. 

Input Keyword giving the number of engineering 
unit specificatins which must immediately 
follow 

Input Keyword giving the number of force 
extra equation specifications which much 
immediately follow 

Input Keyword giving the number of engineering 
unit extra equation specifications which must 
immediately follow 

Uncorrected delta normal force component 
(Main Balance) 

Uncorrected delta normal force component 
(Second Balance) 

Normal force in pounds due to base pressures 
which is to be applied as a correction to 
balance loads 

Balance axis normal force component 

Correct delta normal force component 
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Standard Name 
NFCH 


NFDF 


NFLO 


NFMA 

NFTA 

NFTO 

I.FWA 

NGB 


NGP 


NHVK 


NOOP 


APPENDIX B 

Description 

Normal force in pounds due to chamber 
pressures which is to be applied as a 
correction to balance loads 

Input D-type constant giving normal force 
sting bending deflection constant in 
degree s/pound 

Input Keyword giving the number of flowmeter 
specifications which must immediately 
follow 

Model axis normal force component 

Normal force weight tare component 

Correct total normal force component 

Wind axis normal force 

Input Keyword giving the number of rotation 
specifications from the gravity axis to 
the balance axis which must immediately 
follow. The maximum is 12. 

Input Keyword giving the number of summary 
print specifications which must immediately 
follow. The maximum is 5. 

Input Keyword giving the mamber of reorder 
hookup specifications which must immediately 
follow 

Extra equation type identifier 
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Standard Name 


NPLT 


NPG 


NRTO 


PI 

PGn 


PHIG 

PHII 

PHIS 

PHIW 

PHIZ 


PI 

PIPE 

PM 


Description 

Input Keyword giving the number of normal 
point-by-point plot specifications which 
must immediately follow. The maximum is 4. 
Input Keyword giving the number of 

point-by-point print specifications which 
must immediately follow 
Input Keyword giving the number of pressure 
ratio specifications which must immediately 
follow 

Free stream test section static pressure in PSF 
Input Keyword associated with print 
specifications 

Gravity axis roll angle in degrees 
Initial sting roll deflection angle in 
degrees 

Sting roll deflection angle in degrees 
due to load 
Roll angle 

Roll attitude of the balance with respect 
to gravity when the wind -off zero is 
recorded, in degrees 
Indicated tunnel static pressure in PSF 
Uncorrected static pressure 
Uncorrected delta pitching moment component 
(Main balance) 
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Standard Name 
PM2 


PMB 


PMBA 

PMC 

PMCH 


PMDF 


PMMA 

PMSA 

PMTA 

PMTO 

PMWA 

QINF 

QPR 

Rll 

R21 

• 

R33 


Description 

Uncorrected delta pitching moment component 
(Second balance) 

Pitching moment in inch pounds due to base 
pressures which is to be applied as a 
correction to balance loads 
Balance axis pitching moment component 
Correct delta pitching moment component 
Pitching moment in inch pounds due to 
chamber pressures wliich is to be 
applied as a correction to balance loads 
Input D-type constant giving pitching moment 
sting bending deflection constant in 
degrees/inch-pound 
Model axis pitching moment component 
Stability axis pitching moment 
Pitching moment weight tare component 
Correct total pitching moment component 
Wind axis pitching moment 
Free stream test section dynamic pressure 
in PSF 

Uncorrected dynamic pressure 
The nine elements of the balance 

attitude transformation matrix [r] 
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Standard Name 
REFL 


REYN 

RH 

RHO 

mo 

RM 

RM2 

RMB 


RMBA 

RMC 

•<MCH 


RMDP 


Description 

Input D-type constant giving the length in 
feet for the calculation of Reynolds 
number. If oraitted, the Reynolds number 
per foot is calculated. 

Free stream Reynolds number in millions 
Free Stream relative humidity in percent 

3 

Free stream density in slugs/ft 
Extra equation type identifier 
Uncorrected delta rolling moment component 
(Main balance) 

Uncorrected delta rolling moment component 
(Second balance) 

Rolling moment in inch pounds due to base 
pressures which is to be applied as a 
correction to balance loads 
Balance axis rolling moment component 
Correct delta rolling moment component 
Rolling ^ment in inch pounds due to cheunber 
pressures which is to be applied as a 
correction to balance loads 
Input D-type constant giving rolling moment 
sting bending deflection constant in 
degrees/inch-pound 
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St . d ard Name 
RMMA 

msh 

RUN 

S 

SET 

SETN 

SF 

SF2 

SFB 


SFBA 

SFC 

SFCH 


SFDF 


SFMn 


Description 

Model axis rolling moment component 
Stability axis rolling moment 
Run number 

Inp\it D-type constant giving model reference 
area in square feet 
Extra equation type identifier 
Extra equation type identifier 
Uncorrected delta side force component 
(Main balance) 

Uncorrected delta side force component 
(Second balance) 

Side force in pounds due to base pressures 
which is to be applied as a correction to 
balance loads 

Balance axis side force component 
Correct delta side force womponent 
Side force in pounds due to chamber 
pressures which is to be applied as a 
correction to balance loads 
Input D-type constant giving side force 
sting bending deflection constant in 
degrees/pound 

Flowmeter supercompressibility factor of 
Flowmeter n 
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Standard Name 

Description 

SFMA 

Model axis side force component 

SFSA 

Stability axis side force 

SFTA 

Side force weight tare component 

SFTO 

Correct total side force component 

SFWA 

Wind axis side force 

SIN 

Extra equation type identifier 

SOUT 

Input Keyword specifying the start of output 
SIB data. The data base following 
SOUT will be retained during data base 
collapses. 

SPRn 

Static pressure ratio of Flowmeter n 

SPLT 

Input Keyword giving the number of summary 
plot specifications which must immediately 
follow. The maximum is 4. 

SSC 

Extra equation type identifier 

SUMV 

Extra equation type identifier 

TABl 

Extra equation type identifier 

TAB2 

Extra equation type identifier 

TAB3 

Extra equation type identifier 

TAB4 

Extra equation type identifier 

TAN 

Extra equation type identifier 

TCJ 

Extra equation type identifier 

TCK 

Extra equation type identifier 

TDEW 

Indicated tunnel dewpoint temperature in 

Or 
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Standard Name 

TEST 

TFn 

THEN 

THTn 

TINT 

TT 

VI 

V2 

V3 

VAAV 

VAS 

VASQ 

VAV 

VAVG 

VDS 

VDV 

VEL 

VINE 

Vise 

VMS 

VMMV 

VMV 

VPR 


. ^ 
3 


Description 

Test number 

Flowmeter temperature in Flowmeter n 

Extra equation type identifier 
Normalized temperature of Flowmeter n 
Free stream static temperature in 
Indicated tunnel total temperature in 
Axial force weight tare per unit weight 
Side force weight tare per unit weight 
Normal force weight tare per unit weight 
Extra equation type identifier 
Extra equation type identifier 
Extra equation type identifier 
Extra equation type identifier 
Extra equarion type identifier 
Extra equation type identifier 
Extra equation type identifier 
Extra equation type identifier 
Free stream test section velocity in ft/sec 
Free stream test section absolute viscosity 
in slugs/ft-sec 

Extra equation type identifier 
Extra equation type identifier 
Extra equation type identifier 
Uncorrected velocity 
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Standard 

VSS 

VSSQ 

VSV 

WAF 


WMli 

WM21 

• 

0 

0 

WM33 

WNF 


WPn 

WPNn 

WSF 


WXPM 


Name Description 

Extra equation type identifier 
Extra equation type identifier 
Extra equation type identifier 
Input D-type constant giving model weight 
tare in pounds for balance axial beam. 

A card input overrides the value calculated 
from a tare run. 

The nine elements (WM. J of the model attitude 

ID 

transformation matrix [r 1 

wm 


Input D“type constant giving model weight 
tare in pounds for balance normal beam. 

A card input overrides the value 
calculated from a tare run. 

Weight flow rate of Flowmeter n 
Normalized weight flow rate of Flowmeter n 
Input D-type constant giving model weight 
tare in pounds for balance side beam. A 
card input overrides the value calculated 
from a tare run. 

Input D~type constant giving model x 
moment tare in inch-pounds for balance 
pitch beam. A card input overrides the 


value calculated from a tare run. 
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Standard Name 
WXYM 


WYRM 


WYYM 


WZPM 


WZRM 


XAFB 


Description 

Input D-type constant giving model x 
moment tare in inch-pounds for balance 
yaw beam. A card input overrides the 
value calculated from a tare run. 

Input D-type constant giving model y moment 
tare in inch-pounds for balance roll 
beam. A card input overrides the value 
calculated from a tare run. 

Input D-type constant gxving model y 
moment tare in inch-pounds for balance 
yaw beam. A card input overrides the 
value calculated from a tare run. 

Input D-type constant giving model z 
moment tare in inch-pounds for balance 
pitch beam. A card input overrides the 
value calculated from a tare run. 

Input D-type constant giving model z moment 
tare in inch-pounds for balance roll 
beam. A card input overrides the value 
calculated from a tare run. 

Axial force in pounds due to base pressures 
which is not to be applied as a correction 
to balance loads 
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Standard 

XAFC 

XBAR 

XCAB 

XCAC 

XFLO 

XNFB 

XNFC 

XPMB 


Name Description 

Axial force in pounds due to chamber 

pressures which is not to be applied as a 
correction to balance loads 
Input D-type constant giving x moment transfer 
distance in inches measured in model axis 
system 

Coefficient of axial force due to base 
pressures which contributed to XAFB. 
Coefficient of axial force due to chamber 
pressures which contributed to XAFC 
Input D-type constant giving cross flow 
angle in degrees as a yaw rotation from 
the wind axis to the gravity axis. 

Normal force in pounds due to base pressures 
which is not to be applied as a correci-ion 
to balance loads 

Normal force in pounds due to chamber 

pressures which is not to be applied as a 
correction to balance loads 
Pitching moment in inch pounds due to base 
pressures which is not to be applied as a 


correction to balance loads 
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Standard Name 
XP^C 


XRMB 


XRMC 


XSFB 


XSFC 


XYMB 


XYMC 


YAWG 

YAWI 


Descriptio n 

Pitching moment in inch pounds due to chamber 
pressures which is not to be applied as a 
correction to balance loads 
Rolling moment in inch pounds due to base 
pressures which is not to be applied as a 
correction to balance loads 
Rolling moment in inch pounds due to chamber 
pressures which is not to be applied as a 
correction to balance loads 
Side force in pounds due to base pressures 
which is not to be applied as a correction 
to balance loads 

Side force in pounds due to chamber pressures 
which is not to be applied as a correction 
to balance loads 

Yawing moment in inch pounds due to base 
pressures which is not to be applied as a 
correction to balance loads. 

Yawing moment in inch pounds due to chamber 
pressures which is not to be applied as a 
correction to balance loads 
Gravity axis yaw angle in degrees 
Initial sting yaw deflection angle in 
degrees 
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Standard Name 
YAWS 


YAWW 

YAWZ 


YEAR 


Yn 

YM 

YM2 

YMB 


YMBA 

YMC 

YMCH 


Description 

Sting yaw deflection angle in degrees due to 
load 

Yaw angle in degrees 

Yaw attitude of the balance with respect to 
gravity when the wind-off zero is 
recorded^ in degrees 

Input D-type constant giving y moment 
transfer distance in inches measured in 
model axis system 

Flowmeter expansion factor of Flowmeter n 

Uncorrected delta yawing moment component 
(Main balance) 

Uncorrected delta yawing moment component 
(Second balance) 

Yawing moment in inch pounds due to base 
pressures which is to be applied as a 
correction to balance loads 

Balance axis yawing moment component 

Correct delta yawing moment component 

Yawing moment in inch pounds due to chamber 
pressures which is to be applied as a 


correction to balance loads 
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Standard Name 
YMDF 


YMMA 

YMSA 

YMTA 

YMTO 

YMWA 

ZBAR 


ZERO 

ZMUn 


Description 

Input D-type constant giving yawing mon.ent 
sting bending deflection constant in 
degrees/ inch-pound 
Model axis yawing moment component 
Stability axis yawing moment 
Yawing moment weight tare component 
Correct total yawing moment component 
Wind axis yawing moment 
Input D-type constant giving z moment 
transfer distance in inches measured in 
model axis system 
The constant 0.0 

Flowmeter fluia viscosity of Flowmeter n 
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RTAT CARD INPUT SPECIFICATIONS 

RTAT requires certain input specifications, keywords, and constants 
in order tc work properly. These are usually input in the form of cards 
keypunched from input setup sheets. In certain cases, constants may be 
input through the digital constants panel or through the graphics terminal. 

Two basic types of constants are used by the RTAT task. They are 
classified as SIF-type constants and non-SIF-type constants. SIF-type 
constants are further subdivided into C-type constants and D-type constants. 
When RTAT is processing SIF-type constants, it places the constant name in 

the names disk file and the constant value in the corresponding 

location in the values disk file. The data base will thus 

contain either the name for names records or the value for other types of 

data records. For C-type constants, RTAT searches the names file from the 

beginning for a matching name. If one is found, the corresponding value in the 

values file is replaced with the value from the input card and the search 

in terminated; otherwise the constant is added to the end of the files. 

Since the data base initially consists of all of the input specifications, 

C-type constants must be used with great care. D-type constants are added 

to the end of the files without searching for a matching name. 

When RTAT is processing non-SIF-type constants, it places both the 

constant name and the constant value in both the names and values files. 

Thus, both the name and the value will be contained in the data base for 

both names records and other types of data records. The constant value 

immediately follows the name in the data base. RTAT also uses input keywords 
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and specif ications. A keyword is a non-SIF-type constant name which has a 
definite meaning to RTAT. It may or may not have an associated value. An 
input specification is an ordered sequence of names and/or values with 
definite meaning to RTAT. Specifications may be input in blocks which are 
always identified by a preceeding keyword whose associated value gives the 
number of specifications in the block. A specification block may be 
terminated by a keyword. 

An entire input setup, consisting of constants, keywords, and specifications, 
is collectively called the "input constants" or the "setup." 

The card input specifications of the extra equation capability are 
given separately in APPENDIX D. 

All input cards are punched free field with a blank as a field delimiter. 

An asterisk in the first field signifies that the card contains a comment. 

Names may be optionally bracketed with dollar signs. All names are truncated 
to four characters with the exception of balance names and balance calibration 
dates which are truncated to eight characters. Numeric values may be 
punched as integers or as fixed or floating point constants as desired. 

All numeric values are stored as real numbers. 

In the following card graphics, a vertical line is used to indicate the 


blank or field delimiter. 
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ENGINEERING UNIT EQUATION SPECIFICATIONS 
Engineering unit equations are concerned with converting the raw data 
inputs obtained from OAP into engineering unit quantities. 


^ * CARDS 
NEU_ 


FOR EU SEGfffiNT ] 
f Enter number of EU equations 


OUTPUT INPUT 1 INPUT 2 EQUATION SLOPE INTERCEPT REFERENCE 
NAME NAME NAME TYPE ZERO VALUE VALUE NAME 



, I ; . .1. 1 . 

ENDQ ^ Collapse E.U. constants from ENDQ 


The format of each specification is : 
Field Contents 

1 Output name 

2 Input 1 name 

3 Input 2 name 

4 Equation type and 

zero option 


Description 

The name to assign to the output EU 
value computed by this equation. 

The data base name to be used as the 
primeury input to this equation. 

The data base neune to be used as the 
secondary input to this equation. 

Two concatenated two-character codes 
defining the equation type wd the 
zero option to be used. 
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Field Contents Description 

5 Slope The value of the primary transducer 

sensitivity to be used by this 
equation. 

6 Intercept The value of the offset or secondary 

transducer sensitivity to be used by 
this equation. 

7 Reference name The data base name to be used as a 

reference by this equation. 

Names used to specify Input 1 and Input 2 must exist in the data base 
before the equation using them is executed. They may be defined as output 
names by proceeding engineering unit equaL-lv-^ns, A name used to specify a 
reference value must exist in the data base after all engineering unit 
equations have been executed; hence, it may be defined by a following 
engineering unit equation. 

There are four basic equation types for non-scanivalve data. Each 
equation type may use two wind-off zero options: yes (YE) or no (NO). 

If the zero option code for an equation is YE, the zero value will be 
computed as shown below for each equation type and saved in a special COMMON 
area whenever a wind-off zero record is processed. If the zero option code 
is NO I the zero value is stored as 0.0* 

Equation type code LI designates a linear equation, computed as: 

ZERO » (INPUT1/INPUT2)*SL0PE 

OUTPUT » (INPUT1/INPUT2)*SL0PE+INTERCEPT-ZER0+REFERENCE 
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Equation type code L2 designates a linecu: equation with plus and minus 
slopes, computed as; 

IF (INPUT1/INPUT2) < 0, SENS=» INTERCEPT 
IF (INP0T1/INPUT2) > 0, SENS=SLOPE 
ZERO = (Ii>IPUTl/INPUT2)*SENS 
OUTPUT = (INPUTl/INPUT2)*SENS-ZERO+REFERENCE 
Equation type code AS designates a Kearfott arcsine equation. Equation 
type AS ignores the zero option field and always uses a zero option of YE. 
The equation is computed as; 

ZERO = INPUTl-(lNPUT2/SLOPE)*SIN(INTERCEPTf REFERENCE) 

OUTPUT = ARCSIN((INPUT1-ZEK0)*(SL0PE/INPUT2)) -INTERCEPT 
Equation type code BT designates a linear equation with multiple 
sensitivities. This equation is used with automatic ranging pressure 
transducers such as Baratrons and Baracells. These auto ranging transducers 
may record data in any of seven ranges £Uid change ranges dynamically 
depending on the magnitude of the pressure being measured. Since each range 
has its own calibration slope and intercept, it is necessary to know %diich 
range was used to record a p£u:ticular reading in order to apply the proper 
sensitivity constants. This is accomplished by a second input called the 
range channel. The range channel is specified by the INPUT 2 name while the 
data reading is specified by the INPUT 1 neune. 
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The range channel may be either an analog channel with various millivolt 
levels denoting the ranges, or a digital channel with the actual range 
number available. In either case, the RTAT functi n IRANGE determines the 
range number: 

RANGE = IRANGE (INPUT2) 

where it is required that INPUT2 be available in the data base before the 
BT equation is executed. RANGE must lie in the interval 1 < RANGE <7. If 
the range is invalid, it is assumed to be one. 

The calibration slopes and intercepts are expected to inmediately 
follow the BT equation specification in the data base with the following 
format: 


SLOPEl SL0PE2 . . . SLOPE? Sever; slopes, one per range 

INTERCEPTl INTERCEPT 2 ... INTERCEPT 7 Seven intercepts , one per 

range 


The BT equation therefore is executed as: 

RANGE » IRANGE (INPUT2) 

SLP s SLOPE (RANGE) 

XCPT » INTERCEPT (RANGE) 

ZERO « INPUT1*SLP 

OUTPUT - INPUT*SLP+XCPT-ZERO+REFERENCE 

The zero code for scanivalve data is the port zero option (PO) . 
this may be used with equation types LI, L2, and BT. In the case of 
scanivalves, the hoate port (port zero) reading from each data record is 
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used as the ZERO value for that record. The Output neune must be specified 
as a four-character name with at least the last two ci c^ers numeric. 

The generated names will be sequenced by one for each port. The Input 1 
name is specified as 3n00 for scanivalve n provided Sn is the OAP channel 
name for that analog channel. For the BT equation, a range recording for 
each scanivalve port is required. These range recordings must be contiguous 
and must be in correspondence to the Input 1 port r'^adings. This is 
accomplished ty a specification such as; 

P200 S200 SlOO BTPO 1. 0. PI 

Note that if OAP averaging is utilized, all Baratrons must be set on 
fixed ranges and not permitted to autorange while a test point is being 
recorded . 

Note that an engineering unit equation which does not fit one of the 
above types can probably be handled with the extra equation capeibility invoked 
before force computations. This capability is described in APPENDIX D. 

Number ot Balances Card 

(^BAL I 1 Enter number of balances 

For a force test, this caxd is used to specify the number of balances. 

The current RTAT maximum is two balances, a main balance and a nose balance. 
Note that, due to the design of RTAT, the model attitude computations are 
done within the . . computation modules, therefore, it may be necessary 

to set NBAL to one rather than zero for a non- force test. 
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Balance Interaction Cards 

RTAT permits two methods of entering balance i iteractions . In the 
first method/ the interaction deck is physically included in the setup deck; 
wh' le in the second method, the interaction deck is retrieved from a disk file. 
The card input for the first method is: 

I INTR I 


Index Value Value Value Value Value 



The block of cards following the INTR card is referred to as an I-card 
interaction deck. It is designed to load balance infomation directly into 
a special COMMO .4 block. The Index field specifies the starting location in 

COMMON to load the values contained in the subsequent fields. The I-card 
interaction deck is a standard deck used at Langley vdiich is obtained from 
a group with the responsibility of maintai \ng the balance interaction deck 
library. The contents of the I-card interaction deck are as follows: 

Index Description 

1 Balance identification, treated as a name 

2 Balance calibration date, treated as a name 

3 Balance calibration date in a special integer format 

4 Number of balance components physically defined 

5-10 Array of balance component names. All baiamce matrices 

are arranged in this order. 
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I luit'x 

II 


1.' 
1 ^ 


M 


U> 




«> 


180 - UV> 


18 S-p)l 


VK) 




108-00 1 


00.1-01 S 


OU^ 


ipt ion 

i;poritviiui a nonnal foi co“pi t oin nq moment tv'pe 
l\iliin-'e or a normal force l-nomial foice 2 type balar'-'e 
Oi 0 r of Uilanci' calibration 

Opt ion to translate interactions for initial lotids 
‘Option for one d i scont imious second order interaction 
liuiex in Ku st (’‘ider Inverse times Second Order matrix 
for ill scont inuous second order interaction 
Maximum immbei of iterations allowed for converqence 
InvtMSt' v't NiUma I 1 /.ed Kir:;t Oidt'i' I nt ei\u't ions with main 
diaqonal t'lcnnents of unity (Kirst Order Inverse Matrix) 
Kioduct of Kirsd t>rdei Inverst' Matrix and Nonualii:ed 
Second Ot'iior I nt eract i ons . (First Older Inverse times 
Second Oi di'i Matiix) 

Pojiitive calibration constants 
Ncqat ivt' iMlibration constants 
Percent acciuacy i equirt\l for conviM qonce 
Oalibration prime sensitivities 

Accuracy infoimation used to establish tlie intei'action 
convc'iqtMice cr i ter ia 
Shunt calibration bridqe resistances 

Malance delta weiqht, the balance weiqht not measured 
by the axial lH\\m 
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The card input for the second method is: 
BALANCE NAME CALIBRATION DATE 

r~BAijv r ( 1 


RTAT considers the balance name and calibration date to be eight characcer 
names. RTAT will search a balance history file on the disk for a matching 
balance name and calibration date. If a match if found, RTAT will replace 
the BALA card with an INTR card and a copy of the I-card interaction deck 
obtained from the disk file. No match produces an error message. 

The RTAT balance history file is periodically updated and currently 
contains the following interaction decks: 
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Balance Name 

Calibration Date 

703S 

10/20/77 

718 

08/21/78 

725 

1-7-74 

727 

03/14/78 

728 

04/16/68 

729 

03/30/78 

731 

10/29/65 

731B 

08/01/75 

733 

11/09/77 

73G 

06/13/77 

737-A 

10/20/76 

737/B 

04/21/76 

737/B(A/ 

03/30/76 

738 

09/19/76 

738B 

05/16/73 

739 

06/06/79 

747 

11/07/78 

753 

10/16/78 

804SA 

03/30/79 

804 SB 

03/11/75 

832B 

10/18/77 

832C 

07/17/79 

832D 

10/20/77 

833 

06/27/79 

834 

11/16/77 

835 

11-09-77 

836 

09/20/77 

838 

03/30/78 

839 

07/27/79 

840 

05/29/79 

842A 

05-08-78 

842-B 

11-1-77 

843A 

02/27/79 

843B 

07/05/79 

847 

08/30/76 

847-B 

07/12/77 

UT03-D 

08/12/70 

UT24-50 

05/15/79 

UT24-100 

03/23/72 

UT27-55 

09/23/75 

UT27-55 

09/23/75 

UT27-100 

08/15/75 

UT34A 

3/9/78 

UT34B 

2/14/78 

UT37 

8/3/78 

IR-IO-B 

03/16/77 

IRIO-C 

01/06/72 

1R15 

05/08/72 

IR21 

05/20/77 
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Balance and Model Attitude Cards 

The attitude of the balance and the model (with respect to gravity) are 
specified through an ordered set of Euler rotations. 


*ROTATIONS FROM GRAVITY TO BALANCE 


Enter number of rotations 


NAME AXIS 



*ROTATIONS FROM BALANCE TO MODEL 


Enter number of rotations 


NAME AXIS 
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The format of each specification is: 


Field 

1 

2 


Contents 

Name 

Axis 


Description 

Any valid data base name 
The name PITC^ ROLL, or YAW specifying 
the axis of rotation 


Valid values for NGB and NBM are from zero to twelve 
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Geometrical Constant Cards 



INPUT NAME 

INPUT VALUE 

D 

S 


D 

B 


D 

CBAR 


D 

REFL 


D 

XBAR 


D 

YEAR 


D 

ZBAR 
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Bending Deflection Constant Cards 


INPUT NAME 

INPUT VALUE 

D 

NFDF 


D 

PMDF 


D 

SFDF 


D 

YMDF 


D 

RMDF 
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Blockage and Jet Boundary Correction Cards 


* BLOCKAGE AND 

JET BOUNDARY CORRECTIONS ] 

blk| 

Enter 0 to omit and 1 t 


to apply corrections 



INPUT NAME 

INPUT VALUE 

D 

KWI 


D 

KBI 


D 

J2 


_D_ 

J3 
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Base Pressure Correction Cards 


* BA£ 

jE PRESSURES 

NBAS 

1 Enter 


Enter number of pressures 


OUTPUT 

NAME 


INPUT 

NAME 


AF 


AREAS 

SF 


NF 


RM 


AREAS*ARMS 
PM YM 


FLAG 


The format of each specification is: 


Field 

1 


Contents 
Output name 


Input name 


Area 


AF 


Area 


SF 


Area, 


NF 


Descript l jn 

The ncime to assign to the individual 
output pressure coefficient computed 
from this input pressure 

The data base name of this input 
pressure in pounds/square foot 

The area in square feet to be used in 
computing the axial contribution of 
this pressure 

The area in square feet to be used in 
computing the side contribution of 
this pressure 

The area in square feet to be used in 
computing the normal contribution 
of this pressure 


63 






APPENDIX C 


Field 

6 

7 


8 


Contents 

Area*Ann„ 

RM 


Area*Arm_ 

PM 


Area*Ann„^^ 

YM 


FLAG 


Description 

The area times arm in square feet 

times inches to be uset in computing 
the roll contribution of this pressure 
The area times arm in square feet 

times inches to be used in coii5)uting 
the pitch contribution of this 
pressure 

The area times arm in square feet 

cimes inches to be used in computing 
the yaw contribution of this pressure 
Flag used to control whether all six 
force and moment components are to 
be corrected for the contributions 
defined by this specification card. 

If negative, the components are not 
corrected. 
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Chamber Pressure Correction Cards 


* CHyiBER PRESSUReFI 

N^CBR [ [ Enter number of pressures. 


OUTPUT 

NAME 


INPUT 

NAME 


AF 


AREAS 

SF 


NF 


RM 


AREAS* ARMS 
PM YM 


FLAG 


The format of each specification is : 


Field 


Contents 


Output name 


Input name 


Area 


AF 


Area 


SF 


Area 


NF 


Description 

The name to assign to the individual 
output pressure coefficienv: computed 
from this input pressure 

The data base name of this input 
pressure in PSF 

The area in square feet to be used in 
computing the axial contribution of 
this pressure 

The area in square feet to be used in 
computing the side contribution of 
this pressure 

The area in square feet to be used 
in confuting the normal contribution 
of this pressure 
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Field 

6 

7 


8 


Contents 

Area^Arm 

RM 


Area*Arm_ 

PM 


Area*Arm 

YM 


FLAG 


Description 

The area times arm in square feet 

times inches to be used in computing 
the roll contribution of this pressure 
The area times arm in square feet 

times inches to be used in computing 
the pitch contribution of this 
pressure 

The area times arm in square feet 

times inches to be used in computing 
the yaw contribution of this pressure 
Flag used to control whether all s’x 
force and moment components are to 
be corrected for the contributions 
defined by this specification card. 

If negative, the components are 
not corrected. 
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Pressure Coefficient Array Cards 


* _ PRESSURE CQ^FICIENT ARRA YS 1 

J Enter number of CP arrays • 


OUTPUT INPUT 
NAME NAME NSIZE 





The format of each specification is: 


Field 

1 


Contents 
Output name 


2 Input name 


3 NSIZE 


Description 

The first name to be generated and 
assigned to the first pressure 
coefficient computed from this 
array 

The data base name of the first 
pressure in pounds/square foot in 
this input pressure array 

Th €2 number of pressures in this array 


The second through NSIZE output names will be generated by incrementing 
the trailing digits of the first output name. 

The input pressure array must occupy NSIZE sequential locations in the 
data base. 
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Pressure Coefficient Ratio Cards 


» pressure* COEFFICIENT_MT_IOS' 

' Enter number of Ratio actuations : 

OUTPUT INPUT 

name name scalar name . NSX2E-. 


The format of each specification is : 


Field Contents 

1 Output name 

2 Input name 

3 Scalar name 

4 NSIZE 


Description 

The first name to be generated and 
assigned to the first ratio computed 
from this ari. 

The data base name of the first 
element in the input array 

The data base name of a scalar 
input value 

The number of elements in this array 


The second through NSIZF output names will be generated by incrementing 
the trailing digits of the first output name. 

The input pressure array must occupy NSIZE sequential locations in the 
data base. 
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Flowmeter Cards 


y Flovmoters ' 

LO ' Enter number of flowmeters 

i’!* tyv'^ codes 


‘rUe flowmeter t^ pe codes are as follows ; 


'IVpe Code 


Flow l>\aio fieri h 1 Eiiiiiber 


1 

1* 

s 

0 

7 

8 
9 


ilol 

31d2 

ilo3 

7C31 

7232 

8332 

8331 

8311 

8312 


RTAT contains the calibration information for the nine Venturi- 
type flowmeters available for use in the 7- by 10-foot high speed tunnel 
listed above. 
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Real Time Display Cards 


' * CARDS FOR DISPLAY J 

'xTOSp r" - - 


Enter number ot names to display. 



Collapse display constants from ENDD. 


The format of each specification is; 


Field Contents 

1 Name 

2 Code 


Description 

The data base name of the item to be 
made available for display 
The three-digit thumbwheel number to 
be associated with this data item 


The thumbwheel number associated with a data item must be in the range 
£ code £ 997, a*id must be unique for each item. 




V 
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Point-By-Point Print Cards 


* CARDS FOR PRINT SEGMEHt I 

NPG [ I Eater number of ^specifications 


PG 





■ 


PG 







PG 







PG 







PG 





1 


PG 







PG 







PG 


1 1 





ENDF 


Collappe Print Constants from MDP 


The format of each specification is variable from two to thirteen 
fields: 


Field Contents 

1 Form code 

2 Name 
through 

13 


Description 

A four-character code word describing 
the type of printout desired for 
this specification 
A set of one to twelve data base 
names of data items to be printed 
under this specification 


The form code is a four-character code consisting of the letters "PG” 
followed by two digits. The digits are used to specify the array size 
for data to be printed in columnar format. The digits "00” are used for 
single values. 
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Each specification with form code "PGoo" followed by one to twelve 
names will cause two lines of output to be printed, which contain the 
names and the values corresponding to the specified names. 

If arrays of values (e.g., sets of scanivalve pressures) are to be printed, 
a specification with form code **PGnn** followed by one to six names will 
cause nn lines of output to be printed, each of which will contain the 
values from one location within the requested arrays. Each array 
printed with one specification will be arranged as a column down the page. 

All arrays printed with one specification should contain the Scime number of 
values. The specified array name is considered to be associated with the 
first array location, and the remaining array names are generated by 
incrementing the trailing digits of the specified name. The remaining 
array values are obtained from sequentially following locations. 

Any name on a specification may be followed by a numerical value which 
represents the w.d portion of an F format which will be substituted for the 
default format of F 10.5. 


72 











APPENDIX C 


Summary Print Cards 


^ CARDS FOR RUN SUMMARY SEGMENT 


!NGP 


Enter numter of specifications 


■gp 




— 

— 

[GP . 

GP 

----- 

1 — 

... 

— 

G£_ 



1 

X 


ENDG * Collapse Print Constants from ENDG 


The format of each specification is variable from two to thirteen 
fields: 


Field 

1 


Contents 
Form code 


2 Name 

through 
13 


Description 

A four-character code word describing 
the type of printout desired for this 
specification 

A set of one to twelve data base names 
of data items to be printed under this 
specification 


The form code is a four-character code consisting of the letters "GP” 
followed by two digits. These digits are used to specify how many n 2 unes 
are in that specification. 

The files used to save the summary print data do not permit the saving 
of arrays of data such as pressure coefficient eurrays. 
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Any name on a specification may be followed by a numerical value which 
represents the w.d portion of an F format which will be substituted for the 
default format of F10.5. 
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I nhuk{ I 


Reorder Hookup Cards 
Enter number of names. 



The format of each specification is that each field contains a data 
base name. 

The reoraer cards contain a list of data base names which are relocated 
to a contiguous set of data base locations. These cards are used to rearrange 
a pressure hookup into contiguous locations to set up arrays for plotting 
pressure data. This reordering is done after the data have been printed. 
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Point-By-Point Plot Cards 

All RTAT plotting is based on concurrent plotting of one to four x-y 
plots on a Tektronix 4014 graphics terminal. Each of these plots occupies 
one quadrant of the screen. Each plot is five units by five units. RTAT 
only plots current data, that is, no comparison plotting capability is 
provided. 


r. *P0il^ BY POINT PLOTS 

NPLT 1 Enter number of plots 



[ENDT ] Collapse Plot Specifications 
The format of each specification consists of 10 fields defined as 


follows : 


Field 

1 


Contents 
X Name 


2 X Origin 

3 X Scale 

4 X Label 


5 y Name 


Description 

The data base name of the variable 
for the x-axis 

The value to be assigned to the 
x-^ueis origin 

The scale factor for the x-axis 

A four-character label to use to 
annotate the y-euils 

The data base name of the variable for 
the y-axis 






-i, 


76 



APPENDIX C 




Field 

6 


Contents 
Y Origin 


7 Y Scale 

8 Y Label 


9 


Array size 


10 Symbol code 


Description 

The value to be assigned to the y-axis 
origin 

The scale factor for the y-axis 

A four-character label to use to 
annotate the y-axis 

If greater than 1, both x name and y 
name should contain the same number 
of values. The specified array name 
is considered to be associated with the 
first array location and the remaining 
array values are assumed to sequentially 
follow the first in the data base. 

A code number for the plotting symbol 
to be used chosen from the following 
list: 

Code Symbol 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

The values of NPLT are from zero to four. 
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Circle 

Tilted-plus 

Up-Triangle 

Square 

Star 

Tilted-Square 
Vertical Bar 
Plus 

Up-Arrow 

Down-Arrow 

Dovm-Triangle 
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Additional Point By Point Plot Cards 
If the NPLT specifications were used to plot pressure data, then an 
additional set of specifications may be input to plot additional pressure 
data . 


* ADDITIONAL POINT BY POINT PLOTS 1 

MPLT I [ Enter number of plots 



The format of each specification consists of 10 fields defined as follows: 


Field 

1 


Contents 
X Name 


2 X Origin 

3 X Scale 

4 X Label 


5 Y Name 


6 V Origin 


7 Y Scale 


Description 

The data base name of the variable 
for the x-axis 

The value to be assigned to the 
x-axis origin 

The scale factor for the x-axis 

A four-character label to use to 
annotate the x-axis 

The data base name of the variable 
for the y-axis 

The value to be assigned to the 
y-axis origin 

The scale factor for the y-axis 
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Field 

8 

9 


10 


The values 


Contents 
Y Label 


Array size 


Symbol code 


Description 

A four-chcuracter leUsel to use to 
annotate the y-axis 
If greater than 1, both x name and y 
name should contain the same number 
of values. The specified array 
name is considered to be associated 
with the first array location and iche 
remaining array values are assumed 
to sequentially follow the first in 
the data base. 

A code number for the plotting 
symbol to be used chosen from the 
fi llowing list: 

Code Symbol 


1 

Circle 

2 

Tilted-plus 

3 

Up-Triangle 

4 

Square 

5 

Star 

6 

Tilted-Square 

7 

Vertical Bar 

8 

Plus 

9 

Up-Arrow 

10 

Do«m-Arrow 

11 

Down-Triangle 


of NPLT are from zero to four. 
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Summary Plot Cards 


» SUMMARY p lots"] 

SPLT I [ Enter nun±)er of plots 


I ends] 




Collapse Plot Specifications 


^ j 


I 

I 


The format of each specification consists of 10 fields defined as 
follows; 


Field 


Contents 


X Name 


X Origin 


3 

4 


X Scale 
X Label 


Y Name 


Y Origin 


7 

8 


Y Scale 

Y Label 


Array size 


Description 

The data base name of the variable for 
the x-axis 

The value to be assigned to the x-axis 
origin 

The scale factor for the x-axis 

A four-character label to use to 
annotate the x-axis 

The data base name of the variable 
Jor the y-axis 

The value to be assigned to the 
y-axis origin 

The scale factor for the y-axis 

A four-character label to use to 
annotate the y-axis 

Must be 1 

<4 
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Field Contents Description 

10 Symbol code A code number for the plotting symbol 

to be used chosen from the following 
list: 


Code Symbol 


1 

Circle 

2 

Tiltec-pli’s 

3 

Up-Triangle 

4 

Square 

5 

Star 

6 

Ti 1 ted-Square 

7 

Vertical Bar 

8 

Plus 

9 

Up- Arrow 

10 

Down- Arrow 

11 

Down-Triangle 


The X Name and Y Name names must be availab) -j >. run sunmary print 
files. Si.nce arrays cannot be saved on the rim sunmary files in array from, 
the array size snist not be greater than 1. 

The values of SPLT are from zero to four. 
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RTAT EXTRA EQUATION CARD INPUT SPECIFICATIONS 
Extra Equation Cards 

The extra equation capability may be invoked in two places: at the 

beginning of the NBAL loop before the force computations and at the end 
of the NBAL loop. The former occurs after engineering units and tunnel 
parameters are computed and may be used to compute additional engineering 
unit variables. The latter occurs after UNEEK and may be used to compute 
additional results. The extra equation capability is invoked within the 
NBAL loop. 

The card input for extra equations at the beginning of the NBAL loop 
is* 


toEXU 


Enter number of specification cards 


p 










. .. 









ENDU 

C. llapse specifications 


The contents of the specification cards are described below. 

The card input for extra equations at the end of the NBAL loop is: 


1 - ^ Enter number of specifications cards 


— 








1 


.... 








ENDX ■ 

Collapse specifications 



The contents of the specification ^ards are described below. 
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The extra equation specification cards provide access to a wide range 
of normal arithmetic operations, special algorithm operation:., conditional 
operations, and tabular data operations. Most of the normal and special 
operators will handle either scalar variables or one dimensional arrays. 
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Conditional Operator Clauses 

The extra equations module implements a conditional operator clause 
which consists of an IF expression followed by a THEN expression followed by 
an ELSE expression. An IF expression consists of an IF opt rator and two 
input fields. A THEN expression consists of a THEN operator and any operator 

expression oxco]^t another conditional operator clause. An ELSE expression 
consists of an ELSE operator and any operator expression except another 
conditional operator clause. The IF expression determines, whether t\>^ THEN 
expression or the ELSE expression is the one the user desires to execute. 

Note that, to maintain the integrity of the data base, the non-chosen 
expression is actually executed with the output set to the special value 
-9999.0 and the chosen expression is executed normally. Thus one must use 
conditional operator clauses with care. 
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ABS Operator 

Vector absolute value vector* The specification is: 


OutpuL 0 Input A NSXZE 

Nejne Name 


1 . 


The format of the specification is: 


Field 


Contents 


Description 


Keyword 


Operation code 


Output 0 Name 


The name of the result array which is 


vdiere i goes from 1 to NSIZE. 


added to the end of the data base 


Input A Name 


The data base neune of the input array 


NSIZE 


The number of elements in the o and ^ 


The algorithm executed is 


Oi = la. I 


arrays 
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ACOS Operator 

Vector multiply arccos of scalar. The specification is: 



Output 0 

Input A 

Input B 

NSIZE 


Name 

Name 

Name 


1 ACOS 1 1 


The format of the specification isj 


Field 

1 

2 


3 

4 

5 


Contents 
Keyword 
Output 0 Name 

Input A Name 
Input B Name 
NSIZE 


Description 
Operation code 

The name of the result array \AxLch is added 
to the end of the data base 
The data base name of the input eurray 
The data base name of the input scalar 
The number of elements in the o and a 
arrays 


The algorithm executed is 


o^ * a^ * COS ^(b) 


where i goes from 1 to NSIZE. 
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ASIN Operator 

Vector multiply arcsin of scalar. The specification is: 


Output O Input A Input B NSIZE 
Name Name Name 



The format of the specification is: 


Field 

1 

2 


3 

4 

5 


Contents 
Keyword 
Output O Name 

Input A Name 
Input B Name 
NSIZE 


The algorithm executed is 


o. = a. * SIN"^(b) 

1 X 


Description 
Operation code 

*nie name of the result array which is 
added to the end of the data base 
The data bar > name c2 the input array 
The data base name of the input scalar 
The number of elements in the o wd a 
arrays 


where i goes from 1 to NSIZE. 
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ATAN Operator 

Vector multiply arctan of scalar. The specification is: 

Output o Input A Input B NSIZE 
Name Name Name 



The format of the specification is 


Field Contents 

1 Keyword 

2 Output 0 Name 

3 Input A Name 

4 Input B Name 

5 NSIZE 

The algorithm executed is 


{Description 
Operator code 

The name of the result array which is 
added to the end of the data base 
The data base name of the input 
array 

The data base name of the input scalar 
The number of elements in the o and a 
arrays 




* TAir^(b) 


where i goes from 1 to NSXZE 


ATNZ 


vrtiere i goes from 1 to NSIZE. 
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CAI Operator 

Scalar internal drag axial force coefficient. The specification is? 


Output 0 Input A Input B Input C 
Name Name Name Name 



The foinnat of the specification is; 


Field Contents 

1 Keyword 

2 Output 0 Name 

3 Input A Name 

4 Input B Name 

5 Input C Name 


Description 

Operation code 

The name of the result scalar which is 
added to the end of the data base 

The data base name of the first 
input scalar 

The data base name of the second 
input scalar 

The data base name of the third input 
scalar 


The algorithm executed is 



* {(— ) 
MACH 


* cos(ALPW) * 


1 + 
( 

1 + 


MACH 


2 

c 


( r)} 

MACH 


. b - Pi 
^ QINF 
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where Input A is exit area, Input B is exit static pressure, Input C is 
exit Mach number. S, PI, MACH, QINF, and ALPW are obtained from the 
data base. 

’ inputs must be in compatible units. Note that the result is for a 


single engine. 
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CAS Operator 

Vector multiply cos of scalar add vector multiply sin o£ scale. 
The specification is: 


CAS 


Output O Input A Input B Input C Input D NSIZE 
Name Name Name Name Name 



The format of the specification is: 


Field Contents 

1 Keyword 

2 Output O Neune 

3 Input A Name 

4 Input B Name 

5 Input C Name 

6 Input D Name 

7 NSIZE 


Description 

Operation code 

The name of the result array vAtich is 
added to the end of the data base 

The data base name of the first input 
array 

The data base name of the first input 
scalar 

The data base name of the second 
input array 

The data base name of the second 
input scalar 

The number of elements in the o, and 
c arrays 


i 
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The algorithm executed is 


» a^ * COS(b) + * SIN(d) 

where i goes from ' to NSZZE. 
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CNI Operator 

Scalar internal drag normal force coafficient. The specification is: 


Output Input A Input B Input C 
I'faine Name Name Name 



A 


The format of the specification is: 

Field Contents 

1 Keyword 

2 Output Name 

3 Input A Name 

4 Input B Name 

5 Input C Name 

Ttie algorithm executed i^^ 

o - (») * A . 

's' PI MACH 


Description 

Operation code 

The name of the result scalar vdiich is 
added to the end of the data base 

The data base name of the first input 
scalar 

The data base name of the second 
input scalar 

The data base name of the third input 
scalar 


. 4. £_ 1/2 

^ T 

* SIN(ALPW) * ( — ) 


1 ♦ 


MACH 


i 

■j. 
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where Input A is exit area. Input B Is exit static pressure^ Input .s 
exit Mach number. S, PI, MACH, ind ALPW are obtained frcm the data base. 

The inputs must be in compatible tinits. Note that the result is for a 
single engine. 
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(X)S Operator 

Vector multiply cos of scalar, “Bie specification is: 


fc03 


Output O 
Name 


Input A Input B NSIZE 
Name Name 

I ' ~i I 


The format of the specification is: 


Field 

1 

2 


Contents 
Keyword 
Output O Name 


3 Input A Name 

4 Input B Name 

5 NSIZE 


The algorithm executed is 


Description 
Operation code 

The name of the result 2 urray «diich is 
added to the end of the data base 
The data base name of the input array 
^le data base name of the input scalar 
The number of elements in the o and a 
arrays 


o^ - a^ *COS(b) 

where i goes from 1 to NSIZE. 


PAGE JS 
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CP Operator 

Vector pressure coefficient. The specification is: 


Output 0 Input A Ni^IZE 
Name Name 


CP 



■:] 


The format of the specification is: 


Field Contents 

1 Keyword 

2 Output O Name 

3 Input A Name 

4 NSIZE 


Hie algorithm executed is 


Description 
Operation code 

The name of the result array vAiich is 
added to the end of the data base 
The data base name of the input array 
The number of elements in the o and a 
arrays 


= (a^ - PDQINF 


where i goes from 1 to NSIZE, PI and QINF are obtained from the data 
base. Input A, PI, and QINF must be in compatible units. 
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CSS Operator 

Vector multiply cos of scalar subtract vector multiply sin of scalar. 
The specification is: 


1 


CSS 


Output O Input A Input B Input C Input D NSIZE 
Name Name Name Name Name 



' i 




Olie format of the specification is: 


Fxeld Contents 

1 Keyword 

2 Output O Name 

3 Input A Name 

4 Input B Name 

5 Input C Name 

6 Input D Name 

7 NSIZE 

The algorithm executed is 

« a^^ * COS(b) - * SlN(d) 

«Atere i goes from 1 to NSIZE. 


Description 

Operation code 

The name of the result array vhich is 
added to the end of the data base 

The data base name of the first input 
array 

The data base name of the first input 
scalar 

The data base name of the second 
input array 

The data base name of the second 
input scalcu: 

The number of elements in the o, a, 
and c arrays 
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ELSE Operator 


The specification is: 


I ELSE I 


The format of the specification is: 


Field Contents 


Description 


1 


Keyword 


Operation code 
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EXP Operator 

Vector exponential scalar. The specification is: 


Output O Input A Input B NSIZE 
Name Name Name 



is: 


The format of the specification 

Field Contents 

1 Keyvford 

2 Output 0 Name 

3 Input A Name 

4 Input B Name 

5 NSIZE 

The algorithm executed is 



where i goes from 1 to NSIZE. 


Descritpion 
Operation code 

The name of the result array v4\ich 
is added to the end of the data base 
The data base name of the input array 
The data base name of the input scaleu: 
The number of elements in the o and a 
arrays 
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IFBQ Expression 

If equal expression. The specification is: 


Input A Input B 
Name Name 



The format of the specification is: 


Field Contents 

1 Keyword 

2 Input A Name 

3 Input B Name 

The algorithm executed is: 


Description 
Operation code 

The data base name of the first input 
scalar 

The data base name of the second 
input scalar 


If a » b, execute the THEN expression, 

otherwise execute the ELSE e:^ression. 


or 


A 
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IFGE Expression 

If greater than or equal expression. The specification is: 


Input A Input B 
itjame Nstme 



The format of the specification is: 


Fiel d Contents 

1 Keyword 

2 Input A Name 

3 Input B Name 

Itie algorithm executed is: 


Description 
Operation code 

The data base name of the first input 
scalar 

The data base name of the second 
input scalar 


If a ^ b, execute the TH^ e3q;>ression, 

otherwise execute the ELSE expression. 
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IFGT Expression 

If greater than expression. The specification is: 


Input A Input B 
Name Name 


I IFGT 

I 





The format of the specification is: 


Field Contents 

1 Keyword 

2 Input A Name 

3 Input B Name 

The algorithm executed is: 


Description 
Operation code 

The data base name of the first input 
scalar 

The data base name of the second 
input scalar 


If a > bi execute the THEN expression f 

otherwise execute the ELSE expression. 
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IFLE Expression 

If less than or equal expression. The specification is: 


Input A Input B 
Name Name 


IFLE 


The format of the specification is; 


Field Contents 

1 Keyword 

2 Input A Name 

3 Input B Name 

The algorithm executed is: 


Description 
Operation code 

The data base name of the first input 
scalar 

^e data base name of the second input 
scalar 


If a < b« execute the THEN expression, 

otherwise execute the ELSE expression. 
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IFLT Expression 

If less than expression. The specification is: 

Input A Input B 
Name Name 

lisiL' " 7 ~T I 


The format of the specification is: 


Field 

Contents 

Description 

1 

Keyword 

Operation code 

2 

Input A Name 

The data base name of the first input 



scalar 

3 

Input B Name 

The data base name of the second 



input scalar 


The algorithm executed is: 

If a < b, execute the THEN expression, 

otherwise execute the ELSE eiqpression. 
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IFNE Expression 

If not equal expression. The specification is: 


Input A Input B 
Name Name 



The format of the specification is: 


Field 

Contents 

Description 



1 

Keyword 

Operation code 



2 

Input A Name 

The data base name of 

the 

first input 



scalar 



3 

Input B Name 

The data base name of 

the 

second 



input scalar 




The algorithm executed is: 

If a ^ b, execute the THEN expression^ 

otherwise execute the ELSE e^q^ression. 
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MACH Operator 

Scalar Mach number. The specification is: 



Output 0 Input A 

Input B 


Name Name 

Name 

fuACH 1 


j 


The foncat of the specif icalion is: 

Field Contents 

1 Keyword 

2 Output 0 Nanr 

3 Input A Name 

4 Input B Name 

The algorithm executed is 



Description 
Operation code 

The name of the result scalar idiich is 
added to the end of the data base 
The data base of the first input 
scalar 

The data base name of the second 
input scalar 


where Input A is a static pressure and Input B is a total pressure. 
Input A and Ii^t B must be in compatible units. 
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MFLO Operator 

Scalar mass flow raci.Q« The iJi; 


Output 0 Input A Input B Input C Input D 
Name Name Name Neune Name 



The format of the specification is: 


Field Contents 

1 Keyword 

2 Output 0 Name 

3 Input A NasM 

4 Input B NasM 

5 Input C Name 

€ Input D Nasie 

The algorithm executed is 


Description 
Operation code 

The nasw of the result tlar which is 
added to the end of the data base 
The data base name of the first input 
scalar 

The data base nasM of the second input 

scalar 

The data base nssM of the third input 

scalar 

The data base nasM of the fourth input 
scalar 


■ a * b * c 

® “ RHO * d * VINP 
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where Input A is an exit density, input B is an exit area. Input C is an 
exit velocity. Input D is an inlet area, RHO and VINF are obtained from the 
data base. The inputs must be in compatible units. 
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NOOP Operator 

No Operation. The specification is: 


Pnoop 1 


The format of the specification is: 

Field Jontents Description 

1 Keyword Operation code 

This specification performs no operation. 




I 
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RHO Operator 

Scalar density. The specification is: 


Output o Input A Input B Input C 
Name Name Name Name 



The format of the specif icaion is: 

Field Contents 

1 Keyword 

2 Output O Name 

3 Input A Name 

4 Input B Name 

5 Input C Name 

The algorithm executed is 


Description 

Operation code 

The name of the result scalar vrtiich is 
added to the end of the data base 

The data base name of the first 
input scalar 

The data base name of *'he second 
input scalcur 

The data base name of the third input 
scalar 


5/2 

o = [(a - .379 • PSAT^Pg^)/(53.3 * 32.17 • c)] • [l/(l + — ) ] 





tdiere Input A is a total pressure in psf, Input B is a Mach number f Input C 
is a total temperature in degrees Rankine, and PSAT ^^ is the vapor pressure 
at the dew point temperature TDEW. TDEW is obtained from the data base. 
Output 0 is in slugs/ft^. 

Ill 
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SET Operator 

Reset data base to value. The specification is; 



Input A 
Name 

Value 

nsizE 

[set 







The format of the specification is: 


Field Contents 

1 Keyword 

2 Input A Neune 

3 Value 

4 NSIZE 

The algorithm executed is 

aj^ = value 

where i goes from 1 to NSIZE. Note that a forward search of the data base 
is used to locate Input A Name. This is designed to be executed at the 
beginning of the NBAL loop to permit setting values for any varieUsle which 
the force computations would use a forward seaurch to retrieve. 


Description 
Operation code 

The data base name of em array vdiose 
values are to be reset 
The actual scaleu: value for the new value 
The number of elements in the a array 
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SETN Operator 

Reset data base to name. The specification iss 



INPUT A 
NAME 

INPUT B 
NAME 

NSIZE 

( SETN ' 

1 

L 

_J 





The format of the specification is: 

Field Contents 

1 Keyword 

2 Input A Name 

3 Input B Name 

I NSIZE 

The algorithm executed is 


Description 
Operation code 

The data base name of an array \idiose 
values are to be reset 
The data base name of an array vdiose 
values are to be used as the new 
values 

The number of elements in the a and b 
arrays 




i 

■-! 


where i goes from 1 to NSIZE. Note that a forward search of the data base 
is u«ed to locate Input A N£une and a backward search of the data base is 
jsed to locate Input B Name. This is designed to be used at the beginning 
of the NBAIj loop to permit setting values for any variable which the force 
computations would use a forw^u:d search to retrieve. 
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SIN Operator 

Vector multiply sin of scalar. The specification is: 



Output 0 
Name 

Input A 
Name 

Input B 
Name 

NSIZE 

|SIN 

1 


The format of the specification 

Field Contents 

1 Keyword 

2 Output 0 Name 

3 Input A Name 

4 Input B Name 

5 NSIZE 

The algorithm executed is 

o. = a. * SIN(b) 

1 1 

where i goes from 1 to NSIZE. 


is: 

Description 
Operation code 

The name of the result array which is 
added to the end of the data base 
The data base name of the input array 
The data base name of the input 
scalar 

The number of elements in the o and a 
eurrays 
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SSC Operator 

Vector multiply sin of scalar subtract vector multiply cos of scalar. 
The specification is : 



Output 0 Input A Input B Input C Input D NSIZE 
Nane Name Name Name Name 

. I \ r ~T~ I I 


The format of the specification 

Field Contents 

1 Keyword 

2 Output 0 Name 

3 Input A Neune 

4 Input B Name 

5 Input C Name 

6 Input D Name 

7 NSIZE 


Description 

Operation code 

The name of the result array v4iich is 
added to the end of the data base 

The data base name of the first input 
array 

The data base name of the first input 
scalar 

The data base name of the second 
input array 

The data base name of the second input 
scalar 

The number of elements in the o, a, 
and c arrays 
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The algorithm executed is 


* SIN(b) - * COS(d) 

where i goes from 1 to NSIZE. 
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SUMV Operator 

SuRunation of vector. The specification is; 


Output 0 Input A NSIZE 
Name Name 


SUMV j 


The format of the specification is: 


Field Contents 

1 Keyword 

2 Output 0 Name 

3 Input A Name 

4 NSIZE 


Description 
Operation code 

The name of the result scalar which is 
auided to the end of the data base 
The data base name of the input array 
The number of elements in the array 


The algorithm executed is 


NSIZE 

o * E a 
i«l 
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TABl Operator 

One dimensional table lookup with linear interpolation. The specification 

is: 


Output. 0 
Neime 



Input A NROW 
Name 


The format of the specification is: 


Field 

Contents 

Description 

1 

Keyword 

Operation code 

2 

Output 0 Name 

The name of the result scalar which is 
to be added to the end of the data 
base 

3 

Input A Name 

The name of the input scalar 

4 

NROW 

The number of independent variable A 
values in the following table 

Ihe 

specification of the 

immediately following table ist 


Independent Dependent 
Variable A Variable 

Value Value 






APPENDIX 0 


The format of the specification is: 


Field 

Contents 

Description 

1 

Independent Variable A 
Value 

The value of the independent variable 
corresponding to the Input A neime 

2 

Dependent Variable 
Value 

The value of the dependent vari 2 tble 
corresponding to the Output 0 Name 


The algorithm executed is a one-dimensional table lookup with linear 
interpolation 2 uid extrapolation. 
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TAB2 Operator 

Two dimensional table lookup with linear interpolation. The specification 
is: 


Output Input A Input B NROW NCOL 
Name Name Name 



The fozmat of the specification is: 


Field 

Contents 

Description 

1 

Keyword 

Operation code 

2 

Output 0 Name 

The name of the result scalar «Aiich : 
to be added to the end of the data 
base 

3 

Input A NaM 

The name of the first input scalar 

4 

Input B Name 

The nasm of the second input scalar 

5 

NBOW 

The number of independent variable A 
values in the following table 

6 

NCOL 

The nundMr of independent variable B 
values in the following table 
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The specification of the inmediately following table is; 


Independent Independent Independent 

Variable B Variable B Variable B 

Value VaJue Value 


0.0" ' 


Independent Dependent Dependent Dependent 

Variable A Variable Variable Variable 

Value Value Value Value 



Field Contenta 

1 0.0 

2 NCOL Independent variable B 

values 

The format of the second through 


Description 

Dummy value to force table alignment 
The value of the independent variable 
corresponding to the Input B Name 
NROW rows of the specification isi 


Field Contenta Description 

1 Independent veriebls A The value of the independent variable 

values corresponding to the Input A Name 

2 -^NCOL Dependent variable value The value of the dependent variable 

corresponding to the Output 0 NasM 
The algorithsi executed is a two-dimensional table lookup with linear 
interpolation and extrapolation. 
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TAB3 Operator 


Slope and intercept teUole look’.ip. The specification is: 


Output Input A NROW 
Name Name 


The format of the specification is: 


Contents 


Description 


Key%*ord 


Operation code 


Output 0 Name 


The name of the result scalar which is 


to be added to the end of the data 


Input A Name 


The name of the input scalar 


thm number of independent variable A 


values in the following table 
The specification of the immediately following table is: 


Independent Slope Intercept 

Variable A Value Value 

Value 











✓ 
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The format of the specification is: 


Field Contents 

1 Independent veuriable A 

value 

2 Slope value 

4 Intercept value 


The algorithm executed is 


Description 

The value of the independent variable 
corresponding to the Input A Name 
The value of the slope of a piecewise 
linear fit to the dependent variable 
value corresponding to the Output 0 Ncune 
The value of the intercept of a 

piecewise linear fit to the dependent 
Vcuriable vaiue corresponding to the 
Output 0 Name 


o = a * (Slope) ^ + (Intercept)^ 


where 


(Independent Variable ^ * 5 . (Independent Variable A)^ 
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TAB4 Operator 


One-dimensional ted>le lookup with exact match. 


The specification is: 



Output Input A NROW 
Name Name 



The format of the specification is: 


Field 

Contents 

Description 

1 

Keyword 

Operation code 

2 

Output 0 Name 

The name of the result scalar which 
to be added to the end of the data 
base 

3 

Input A N«une 

The name of the input scalar 

4 

NROW 

The number of independent variable A 
values in the following table 


The specification of the immediately following table is: 


Independent 
Variable A 
Value 


Dependent 

Variable 

Value 
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The format of the specification is: 


Field 


Contents 


Description 


1 Independent variable A The value of the independent variable 

corresponding to the Input A name 

2 Dependent variable value The value of the dependent variable 

corresponding to the Output O Name 

The algorithm executed is a one-dimensional table lookup with one to 
one correspondence. If no match is found, an eiror message is generated. 
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TAN Operator 

Vector multiply tan of scalar. The specification is: 


Output 0 
Name 



Input A Input B NSIZE 
Name Name 


The format of the specification is: 


Field Contents 

1 Keyword 

2 Output O Name 

3 Input A Name 

4 Input B Name 

5 NSIZE 

The algorithm executed is 


Description 
Operation code 

The name of the result uray which is 
added to the end of the data base 
The data base name of the input 
array 

The data base name of the input 
scalar 

The number of elements in the o and a 
arrays 


Oj^ ■ a^ * TAN(b) 
vAiere i goes frcmi 1 to NSIZE. 
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TCJ Operator 

Vector Iron/Constantan thermocouple temperature- The specificati.^ji xs: 


OUTPUT 0 INPUT A 

NAM E NA ME N SIZE 

1‘TCJ .I __ ,L ...I, ZU 


The format of the specification is: 


Field Contents 

1 Keyword 

2 Output 0 Name 

3 Input A Name 

4 NSIZE 


Description 
Operation code 

The name of the result array vAiich is 
added to the end of the data base 
The data base name of the input 
2 urray 

The number of elements in the o and a 
^u:rays 


The algorithm executed is a linear interpolation in an iron/constantan or 
iron/copper-nickel thermocouple table (a Type J thermocouple) for a cold 
junction at 32^-. The Input A is in absolute millivolts and the Output 0 
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TCK Operator 

Vector Chromel/Alumel thermocouple temperature. The specification is: 


OUTPUT O 
NAME 


INPUT A 
NAME 


NSIZE 


The format of the specification iss 


Field 


Contents 


Description 


Keyword 


Operation code 


Output 0 Name 


The name of the result array which is 


added to the end of the data base 


Input A Name 


The data base name of the input array 


NSIZE 


The number of elements in the o and a 


arrays 

The algorithm executed is a linear interpolation in a chromel/alumel or 
nickel-chromium/nickel-aluminum thermocouple tcdile (a Type K thermocoiqple) 
for a cold junction at 32°F. The Input A is in absolute millivolts and the 
Output O is in ®F. 


’* . 0 
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THEN Operator 


The specification is: 


I THEN I 


The format of the specification is; 


Field Contents 


Description 


1 


Keyword 


Operation code 
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VAAV Operator 

Vector add vector add vector. The specification is: 


VAAV 


Output 0 
Name 



Input A Input B Input C NSIZE 
Name Name Name 


The format of the specification is: 


Field Contents 

1 Keyword 

2 Output 0 Name 

3 Input A Name 

4 Input B Name 

5 Input C Name 

6 NSIZE 

The algorithm executed is 


Description 

Operation code 

The name of the result array tdiich is 
added to the end of the data base 

The data base name of the first input 
array 

The data base name of the second 
input array 

The data base name of the third input 
array 

The number of elements in the o, a, b, 
and c arrays 


•i * *>1 * °i 


where i goes from 1 to NSIZE. 
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VAS Operator 

Vector add scalar. The specification is: 


Output 0 Input A Input B NSIZE 
Name Name Name 



The format of the specification is: 


Field Contents 

1 Keyword 

2 Output 0 Name 

3 Input A Name 

4 Input B Name 

5 NSIZE 


The algorithm executed is 


Description 

Operation code 

The name of the result array which is 
added to the end of the data base 

The data base name of the input 
array 

The data base name of the input 
scalar 

The number of elements in the o and a 
arrays 


°i 


+ b 


where i goes from 1 to NSIZE. 
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VASQ Operator 

Sum of squared vectors. The sp:icification is: 


Output 0 Input A Input B NSIZE 
Name Name Name 

I VASQ y I 


Die format of the specification is: 


Field Contents 

1 Keyword 

2 Output 0 Name 

3 Input A Name 

4 Input B Name 

5 NSIZE 

Die algorithm executed Is 


Description 

Operation code 

The name of the result array which 
is added to the end of the data 
base 

The data base name of the first input 
array 

The data base name of the second 
input array 

The number of elements in the o, a» 
and b arrays 


°i 


2 

“i 



where i goes from 1 to NSIZE. 
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VAV operator 

Vector add vector- The specification is: 


Output 0 Input A Input B 
Name Name Name 


NSIZE 


L-_. 


Field 


Contents 


Description 


Keyword 


Operation code 


Output 0 Name 


The name of the result array which is 


added to the end of the data base 


Input A Name 


The data base name of the first 


input array 


Input B Name 


The data base name of the second 


input array 


NSIZE 


The number of elements in the o. 


a« and b arrays 


The algorithm executed is: 


Oi - a^ + bi 


where i goes from 1 to NSIZE. 
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VAVG Operator 

Vector average. The specification is: 


Output 0 Input A Input B NSIZE 
Name Name Name 

[vAVG I I 


The format of the specification is: 


Field Contents 

1 Keyword 

2 Output 0 Name 

3 Input A Name 

4 Input B Name 

5 NSIZE 

The algorithm executed is 


Description 

Operation code 

The name of the result array which 
is added to the end of the dati 
base 

The data base name of the first input 
array 

The data base name of the second 
Input array 

The number of elements in the o, a, 
and b arrays 


o^ - (a^ + b^)/2 
sdiere i goes from 1 to NSIZE. 
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VDS Operator 

Vector divide scalar. The specification is: 


Output 0 Input A 
Name Name 


Input B NSIZE 
Name 


[vDS I I 



The format of the specification 

Field Contents 

1 Keyword 

2 Output 0 Name 

3 Input A Name 

4 Input li Name 

5 NSIZE 

The algorithm executed is 


Is: 

Description 
Operation code 

The name of the result array which is 
added to the end of the data base. 
The data base name of the input 
array 

The data base name of the input 
scalar 

The number of elements in the o and a 

arrays 


tdiere i goes from 1 to NSIZE. 
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VDV Operator 

Vector divide vector. The specification is: 


[vDV 


Output 0 
Name 


Input A Input B NSIZE 
Name Name 


The format of the specification is: 


Field Contents 

1 Keyword 

2 Output O Name 

3 Input A Name 

4 Input B Name 

5 NSIZE 

The algorithm executed is 


Description 

Operation code 

The name of the result array which is 
added to the end of the data base 

The data base name of the first 
input array 

The data base name of the second 
input array 

The number of elements in tlie o, a, 
and b arrays 


Oj . .^/bj 


tihere i goes from 1 to N81ZB 
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VEL Operator 

Scalar speed. The specification is: 


Output 0 Input A Input B 
Name Name Name 


“I 


The format of the specification is: 


Field Contents 

1 Keyword 

2 Output 0 Name 

3 Input A Name 

4 Input B Name 

The algorithm executed is 

o = jl.h * 53.3 * 32.1? • 


Description 
Operation code 

The name of the result scalar which is 
added to the end of the data base 
The data base name of the first input 
scalar 

The data base name of the second 
input scalar 


2 

* ( 1/(1 + 1 -)) 


where Input A is a Mach number and Input B is a total temperature in degrees 
Bankine. Output 0 is in feet per second. 
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VMMV Operator 

Vector multiply vector multiply vector# The specification is: 


Output o Input A Input B Input C NSIZE 
Name Name Name Name 

jVMMV I ’ [ 



The format of the specirlcation is: 


Field Contents 

1 Keyword 

2 Output 0 Name 

3 Input A Name 

4 Input B Name 

5 Input C Name 

6 NSIZE 


The algorithm executed is 


Description 

Operation code 

The name of the result array «Aiich is 
added to the end of the data base 

The data base name of the first 
input array 

The data base name of the second 
input array 

The data base n 2 une of the third 
input array 

The number of elements in the o, 
a, b, and c arrays 


A * h ^ c 
*i **i ®i 


where 1 goes from 1 to NSIZE. 


138 


APPENDIX D 


VMS Op. :tor 

Vector multiply scalar. The specification is: 



Output 0 

Input A 

Input B NSIZE 

Name 

Name 

Name 

i 


1 \ '_1 ] 


The format of the specification is: 

Field Contents 

1 Keyword 

2 Output 0 Name 

3 Input A Name 

4 Input B Name 

5 NSIZE 

The algorithm executed is 

o. * a. * b 
1 1 

where i goes from 1 to NSIZE. 


Description 

Operation code 

The name of the result array v;hich is 
added to the end of the data base 

The data base name of the input 
array 

The data base name of the input 
scalar 

The momber of elements in the o and a 
arrays 
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t 

i 

I VMV Operator 


Vector multiply 

vector . 

The specification is: 

Output 0 

Input A 

Input B NSIZE 

Name 

Name 

Name 

[vMV 1 1 


1 1' .1 


The format of the specification is: 


Field Contents 

1 Keyword 

2 Output 0 Name 

3 Input A Name 

4 Input B Name 

5 NSIZE 

The algorithm executed is 


Description 

Operation code 

The name of the result array which is 
added to the end of the data base 

The data base name of the first 
inpuc array 

The data base name of the second 
input curray 

The number of elements in the o, 
a, and b arrays 



where i goes from 1 to NSIZE. 
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VSS Operator 

Vector siibtract scalar. The specification is: 


Output 0 Input A Input B NSIZE 
Name Name Name 




The format of the specification is: 


Field Contents 

1 Keyword 

2 Output O Ncune 

3 Input A Name 

4 Input B Name 

5 NSIZE 

The algorithm executed is 


Description 

Operation code 

The ncune of the result array which is 
added to the end of the data base 

The data base name of the input 
eurray 

The data base name of the input 
scalar 

The number of elements in the o and a 
arrays 


°i ■ *i - “ 


where i goes from 1 to NSIZE. 


A 


V \ 

Y 


/; 


- ' , 
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VSSQ Operator 


Difference of squared vectors. 

Output 0 Input A 
Name Name 

pSSQ 


The format of the specification is: 

Field Contents 

1 Keyword 

2 Output 0 Name 

3 Input A Name 

4 Input B Name 

5 NSIZE 

The algorithm executed is 



where i goes from 1 to NSIZE. 


The specification is: 

Input B NSIZE 
Name 

] 


Description 

Operation code 

The name of the result array viiich is 
added to the end of the data base 

The data base name of the first input 
array 

The data base name of the second 
input array 

The number of elements in the o, a, 
and b arrays 
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VSV Operator 

Vector subtract vector. The specification is: 

Output 0 Input A Input B NSIZE 
Name Name Name 

j VSV I 



The format of the specification is: 


Field Contents 

1 Keyword 

2 Output 0 Name 

3 Input A Name 

4 Input B Name 

5 NSIZE 

The algorithm executed is 


Description 

Operation code 

The name of the result array v^ich is 
added to the end of the data base 

The data base name of the first 
input array 

The data base name of the second 
input array 

The number of elements in the o, a., 
and b arrays 


o 


i 



where i goes from 1 wO NSIZE. 
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RTAT INTERACTIVE CALIBRATION 
SPECIFICATIONS AND EXAMPLES 

Once the calibration data are acquired, RTAT contains the capeibility to 
work up the following types of calibrations: 

(a) Lineeu: transducers 

(b) Kearfott accelerometers 

(c) Balance spans 

(d) Sting bendings 

The interactive device- may be either the teletype or the Tektronix 4014 
graphics terminal. The Tektronix should be set in small character size with 
AUTO PRINT on and MARGIN CONTROL 1 on. The Tektronix hard copy unit should be 
on. 

The general steps in a simple calibration session are: begin calibration 

session, specify type of calibration and options, specify PROCESS directive, 
specify VOID name directive, specify NAME channel directive, specify VOID 
point directive, specify CHANGE data directive, conputer calibration, emd 
output results. The formats of the directives are given below. 
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Calibration Type Directive 

The specification is: 


TYPE OPTION 1 OPTION 2 OPTION 3 

I I I 1 


The format of the specification is: 


Field Contents Description 

1 Keyword XDCR indicates linear transducer 

calibration. KFTT indicates 
Kearfott arcsine calibration. 

BSPN indicates balance span 
calibration. BNDG indicates sting 
bending calibration. 

2 Option 1 flag A zero or blank turns the option off. 

A one turns the option on. If on, 
the user must identify the dependent 
and independent veu:iad>les. If 
off, the dependent variable is the 
first neune in the NAME channel 
directive and the following names 
are the independent variables in the 
order in which the calibrations 
are to be performed. 
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Option 2 flag 


Option 3 flag 


Field Contents Description 

A zero or blank turns the option off. 

A one turns the option o:i. If 
on, all communication will be via 
the card reader (not recomnended) . 

If off, coinnunication will be via 
the interactive device. 

A zero or blamk turns the option off. 

A one turns the option on. If on, 
the calibration will be worked up 
before corrections are requested. 

If off, corrections will be solicited 
before the first workup. 

Option 1 may not be turned on for types B5 :pn and BNDG. Option 3 may not be 
turned on for type BSPN. I4ie usual selection is all three options turned 
off, thus normally only the first field needs to be specified. 
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K 


I 


Process Directive 


The specification is: 


NAME 1 VALUE 1 NAME 2 VALUE 1 NAME 1 VALUE 2 NAME 2 VALUE 2 

1 p«gg-i - 1 I 1 1 thro) 1 I 1 


The format of the specification is: 


Field 

Contents 

Description 

1 

Keyword 

Directive name 

2 

Name 1 

First calibration file vjuricdble name 

3 

Value 1 

Starting value of first calibration 
file variable name 

4 

Name 2 

Second calibration file variable name 

5 

Value 1 

Starting value of second calibration 
file varlad>le name 

6 

THRU 

Keyword, inclusive through 

7 

Name 1 

First calibration file variable mane 

8 

Value2 

Ending value of first calibration 
file variable name 

9 

Name 2 

Second calibration file variable naaie 

10 

Value 2 

Ending value of second calibration 
file variable name 
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Name 2 Value 1 and Name 2 Value 2 are optional on tdxis specification. 

Name 1 Value 1 THRU Name 1 Value 2 are also optional on this specification. 
This directive is used to select that portion of the calibration input file 
which the calibration session will use. The basic command requires only 
the first fif Id and uses the entire calibration input file. 
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VOID Name Directive 


The specification is: 


NAME VALUE VALUE VALUE VALUE 

I VOID 1 I I THRU I I I THRU I I “1 1 


The fonnat of the specification is: 


Field 

Contents 

Description 

1 

Keyword 

Directive name 

2 

Name 

Calibration file variab .e name 

3 

Value 

Starting value of calibration file 
variable name 

4 

THRU 

Keyword, inclusive through 

5 

Value 

Ending value of calibration file 
variable name 

6 

Value 

Starting value of calibration file 
variable name 

7 

THRU 

Keyword, inclusive through 

8 

value 

Ending value of calibration file 
variable name 

9 

-1 

Directive terminator 


THRU is optional in this specification. Value THRU Value uy be repeated 
as often as required in this specification. The null form of this directive 
has two forms: 
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(a) the directive rame and directive terminator, and 

(b) the directive terminator. 

This directive voids for the entire calibration session. 
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NAME Chemnel Directive 

The specification is: 

NAME 1 NAME 2 TPl TP2 

1 Nj^ 1 ) THRU i I 1 t -1 y 


The format of the specification is: 


Field 

1 

2 

3 

4 

5 

6 
7 


Contents 
Keyword 
Name 1 
THRU 
Name 2 
TPl 
TP2 
-1 


Description 
Directive name 

First calibration file variable name 

Keyword, inclusive through 

Second calibration file variable name 

Starting Test Point number 

Ending Test Point number 

Directive terminator 


TPl and TP2 are optional in this specification. THRU is optional in this 
specification. Name 1 THRU Name 2 may be repeated as often as required in chis 
specification. NAME is also optional in this specification. This directive, 
in conjunction with the option 1 flag of the calibration type directive, 
specifies the names of the dependent and independent variables to be read 
from the calibration input file for each test point. 
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VOID Point Directive 

The specification is: 


NAME VALUE 1 VALUE 2 

I VOID { I i THRU I 1 -1 1 


The format of the specification is: 


Field "ontents 

1 Keyword 

2 Name 

3 Value 1 

4 THRU 

5 Value 2 

6 -1 


Description 
Directive name 

Calibration file variable name, for 
exan^le, TP for test point number 
Starting value of calibration file 
variable name 

Keyword r inclusive through 
Ending value of calibration file 
variable name 
Directive terminator 


THRU is optional in this specification. Value 1 THRU Value 2 may be 
repeated as often as required in this specification. 

The null form of this directive has three forms: 

(a) the directive name and directive terminator; 

(b) the directive terminator; and 

(c) the directive name. 

During a calibration session, the void points accumulate. The input of the 
directive name alone will clear out all previous void points. 
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CHANGE Data Directive 

The specification is: 


NAME TP VALUE 

I CHAN 1 

L 1 

L -1 J 

The 

format of the specification 

is: 

Field 

Contents 

Description 

1 

Keyword 

Directive name 

2 

Name 

Calibration file variable 

3 

TP 

Test Point number 

4 

Value 

New value to be used 

5 

-1 

Directive terminator 


Name TP Value may be repeated as often as required in this specification* 

The null form of this directive has three forms: 

(a) the directive name and directive terminator? 

(b) the directive terminator; and 

(c) the directive name. 

During a calibration session, the change points accumulate. The input of the 
directive name alone will clear out all previous change points. 
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The specification is: 


RESTART Directive 


I REST t 


The format of the specification is: 

Field Contents Description 

1 Keyword Directive name 

This directive can be input ciny time a directive is requested. This directive 


causes RTAT to re-ini tialize the calibration session. 
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Delete Input Line Directive 
The specification is: 

BLANK 

1 ■ I. , D i 

The format of the specification is: 

Field Contents Description 

1 Blank Blank 

2 Keyword Directive Name, the letter D 

This directive may be appended to any input line. This directive causes 
the current input line to be deleted and permits the user to re-enter that 
input line. 
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Linear Transducer Calibration (XDCR) 

The linear transducer calibration type is used with absolute and 
differential linear transducers such as Statham pressure gages. 

Multiple transducer calibrations may be performed during a calibration 
session. The calibration input data must have been previously recorded and 
saved on the RAD file CALMV. The load channel is usually defined on the 
Digital Constants Panel and is therefore also recorded with the transducer 
readings on CALMV. The load values should be input in the correct units. 

For the linear transducer calibration, the object is to determine the 
SLOPE and INTERCEPT which appear in the following equation: 

Dependent variable = SLOPE * Independent variable + INTERCEPT 

where a least square fit is performed using all data points. 

The output of the calibration consists of a listing of the input data, 
the calibration constants, the deviation of each point from the calibration 
equation, the root mean square error, the point number of the point having 
maximum deviation, and the average power supply voltage. 
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Kearfott Accelerometer Calibration (KPTT) 

The Kearfott accelerometer calibration is designed specifically for the 
special arcsine equation used for Kearfott accelerometer attitude sensors. 

Multiple Kearfott calibrations may be performed during a calibration 
session. The calibration input data must have been previously recorded and 
saved on the RAD file CALMV. The load channel is usually defined on the 
Digital Constants Panel amd is therefore also recorded with the accelerometer 
readings on CAI/IV. The load values should be input in the correct units. 

For the Kearfott calibration, the object is to determine A, C, and 
PHIO which appear in the following equation: 

Dependent variable = arcsin - PHIO 

where a least square fit is performed using all data points. 

The output of the calibration is a listing of the input data, the 
calibration constants, the deviation of each point from the calibration 
equation, the root mean square error, the point number of the point having 
maximum deviation, and the average power supply voltage. 
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Balance Span Calibration (BSPN) 

The balance span calibration is designed specifically to work up 
balance shunt spans and compute balance sensitivity constants* 

Only one balance span may be performed during a calibration session. 

The calibration input data must have been previously recorded and saved on 
the RAD file CALMV. The span codes are usually defined in a digital channel 
set up for thumbwheel input and are therefore also recorded on CALMV. 

The balance span calibration requires the six balance components in 
consecutive analog channels in the order NF,AF^PM,RM,YM,SF, the power supply 
voltage channel, and the span code channel. 

The span code word identifies the analog channel number for NF, 
identifies the balance component being spanned, and identifies the type of 
span load* The span code word is considered an integer as follows: 


Span Code Word Value 

units digit 1 

2 

3 

tens digit 1 

2 

3 

4 

5 

6 

thousands digit and 
hundreds digit combined 


Description 
zero load 
span load 
return zero load 
normal span 
axial span 
pitch span 
roll span 
yaw span 
side span 

analog channel number of 
normal force channel 
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If any portion of the balance span is unsatisfactory, a complete new 
balance span must be recorded. 

If the balance span is satisfactory, the user has the option to apply the 
span to the balance laboratory calibration to obtain the sensitivity 
constants as installed in the tunnel. To do this, the user must know the 
current values of the la±)oratory prime sensitivity constants, the laboratory 
span, and the ledx>ratory power supply voltage. For convenience, these are 
usually punched on cards in the order noted and in the expected balance 
component order. 

RTAT then calculates the tunnel sensitivity constants from: 

Tunnel sensitivity constant = (Laboratory sensitivity constant) 

★ / Laboratory span\ 

I Tunnel span / 


where 


^ _ 1 /Tu nnel Span Load Reading 

ne span - Power Supply Voltage 


1 J^unnel Zero Load Reading 

2 Tunnel Power Supply Voltage 


Tunnel Return Zero Load Reading 
Tunnel Power Supply Voltage 


(Laboratory Power Supply Voltage) 
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The output of the celibration eonaists of a listing of the input data 

the tunnel spans, the tunnel sensitivity constants, and the average tunnel 
balance voltage , 
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Sting Bending Calibration (BNDG) 

The sting bending calibration is designed to work up combined sting 
emd balance deflection constants for either single components or paired 
components . 

Multiple bending calibrations may be performed during a c^Alibration 
session. The calibration input data are normally previously recorded 
manually and punched on cards for input to RTAT. 

The format of the calibration input deck is a header card followed by 
data cards followed by one or more trailer cards. 

The specification of the header card is : 

Btype 


The format of the specification is: 

Field Contents Description 

1 Btype Keyword 

NF indicates normal force bending 
PM indicates pitching moment bending 
SF indicates side force bending 
YM indicates yawing moment bending 
RM indicates rolling moment bending 
LONG indicates combined normal force and 
pitching mcHnent bending 
LAT indicates combined side force and 
yawing moment bending 
ROLL indicates rolling moment bending 
161 
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Note that the names of the independent variables are determined from the 
Btype. 

The specification of the bending data cards is: 


TP LOCATION LOAD DEG MIN SEC 

I — I I r :' : T" : T : i 


The format of the specification is; 


Field 

1 


Contents 

TP 


Description 

An ascending point number which will 
be used by PROCESS. VOID, and CHANGE 
DIRECTIVES 


2 LOCATION 


3 LOAD 

4 DEG 

5 MIN 

6 SEC 


The location of the applied load with 
respect to the bal 2 mce moment 
center in inches 
The applied load in pounds 
The degrees portion of the angle 
The minutes portion of the angle 
The seconds portion of the angle 


Note that the bending calibration uses the nasie DEFL for the 
variable which it computes as follows; 


dependent 


DEFL • DEG + 


MIH 

60 


SEC 

3600 
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The bending calibration is not aware of trigonometric indentities such 
as 360° ■= 0°. (The bending calibration is actually a special type of 
linear transducer calibration.) Care must be exercised in specifying the 
deflection angles in order to obtain correct results. In view of the 
definition of DEFL, an angle such as 358°20'50" should be input as the DEG 
MIN SEC string -1 -39 -10. 

The specification of the trailer card is s 


I E ND I 


The format of the specification is: 

Field Con tents Description 

1 Keyword One END card terminates a Btype . f a'. 

END cards terminate all bending input. 

For the bending calibration, the object is to determine the values of 
CO, Cl, and C2 which appear in the following equation: 

Dependent variable = CO + Cl * (Independent variable 1) 

C2 * (Independent variable 2) 

where a least squeure fit is performed over all data points. 

The output of the calibration consists of a listing of the input data, 
the calibration constants, the deviation of each point from the calibration 
equation, the root mean square error, emd the point number of the point 
having maximum deviation. 
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Sample Interactive Calibration Sessions 
The following examples represent simple Interactive calibration sessions 
using the Tektronics graphics terminal. The stylized arrow Is the system 
prompt character fur the graphics terminal. The user supplied Input follows 
the prompt character emd is terminated by a carriage return. 

The following is an Interactive linear transducer calibration session. 



IHPuT CAuIhRATION type ( »xnCR» » *KFTr» . »d; 40 G» ♦ ••^SPN* ) AND OPTIONS 
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input change info 

channel '*AME*P0INT NO.tNE'K VALUE 


cUWkEnTlY using aDCK channels LOAU AND S30U 
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The following is an interactive balance span session. 
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The following is an interactive bending calibration session 
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APPENDIX F 


RTAT SAMPLE SETUP DECKS 

The following examples represent simple RTAT setup decks for several 
types of tests. 

These examples are for the tunnel setup only. The setup for the static 
room would consist of a set of input cards similar to the ones following 
the RBM job control cards but specifying the appropriate static room 
channels. Note that OAP can be setup for the tunnel and static room in 
any order. Similarly RTAT can be setup for the tunnel an static room in 
any order. Care should be exercised to avoid conflicts between the tunnel 
and the static room RTAT setups in the ucse of Tektronix 4014 graphics 
terminal. Note that all display thumbwheel assignments must be unique. 

The following is an RTAT setup deck for a pressure model test. 
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^ The following is an RTAT setup deck for a single balance force model test. 
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The following is an RTAT setup deck for a dual balance force model 
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APPENDIX G 


RTAT COMPUTATIONAL SPECIFICATIONS 
RTAT Raw Data Input 

Raw data values are obtained from the data acquisition system by the OAP, 
which averages them over a test point and places the results in an averaged 
record buffer. RTAT obtains these values and places them in the dita base. 

The following conventions are used in assigning names to the raw data 
values: 

(a) Digital channel data are designated by the letter "D" followed by 
the channel number; e.g., D151 and D166. 

(b) Analog channel data (non-scanivalve data only) are designated 
by the letter "A" followed by the channel number; e.g., A1 ^md A47. 

(c) Scanivalve data are designated by the letter "S" followed by a 
three-digit number giving the valve number and port number; e.g., S237 for 
valve 2, port 37. The first two characters, e.g., S2, are obtained from 
the OAP channel neune t 2 d}le. 

(d) Tachometer channel data are designated by the letter "T" 
followed by the chemnel number; e.g., T211. 

(e) The letter "V” followed by a three digit channel number has been 
reserved for the future for a Vidar data chemnel. 

(f) The first four characters of the Digital Constant Panel names 
defined in the OAP setup cards. 

For analog data (including scanivalve data) , RTAT converts the raw data 
to uncorrected millivolts and then applies system calibrates to obtain 
corrected millivolts which it places in the data base. 
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Engineering Units 

RTAT expects the Keyword NEU to be availed>le in the data base. If 
NEU is not found, or if it is less than 1, no engineering unit conversion 
equations will be executed. 

The engineering unit conversion equations will handle 2d3solute and 
differential transducers having linear calibration equations or linear 
equations with plus and minus slopes. They will handle Kearfott arcsine 
equations and Beuratron multiple ramge equations. They will also handle 
scanivalve multiplexed differential pressure transducers. 

The engineering unit input specifications are given in APPENDIX C with 
additional details on the equations inplemented. 

The results obtained through executing the engineering unit equations 
are added to the end of the data base. Then the engineering unit specifications 
are collapsed out of the data base. 

Tunnel Parameters 

RTAT expects to find the following data items in the data base: 

PI, HI, TT, TDEW, and REPL. 

The following equations are used to compute the tunnel flow parameters. 

S<xne of the numeric constants are derived from tunnel calibrations and apply 
only to the Langley 7- by 10-foot high speed tunnel. (See references 1 and 2.) 
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The free stream test section static pressure is computed as: 


For (^) < 0.62474 

HX 




PI = [1.051280* (Ji) - 0.030756]* HI 


For 


A 

'hi 


> 0.62474, 


PI = [0.995865*(^) + 0.004020]* HI 


The relative humidity is confuted as 


VAPl = PSAT(TDEW) 


VAP2 « PSAT(TT) 




where PSAT is a function vfhich gives the saturation vapor pressure of water 
in air at a given tonperature. 
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P'ree stream test section Mach number is computed as : 


MACH = y5*[(Hl/Pl)^^^ - 1] 


Free stream test section dynamic pressure is computed as: 


QIKF = 0.7*P1*MACH^ 


Free stream test section density is computed as: 


RlIO 


2.5 


1 


1 + 0.2 « MACH^ 

*(HI - .379 « VAPl) 


17li*.87i»2»TT 


P’ree strerun test section static temperature is computed as: 


TIhT = TT/(1 0.2 « MACH^) 


Reynolds number in millions based on the length REFL is computed as: 


VINP » ^ * QINF/RHO 


Vise - 2-2T * (TINF ^59.688)^^^ * lO:^. 
TINF + 658.286 


REYN 


s RHO * VINF * REFL , 


Vise 


Note that the default value for REEL is unity. 
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RTAT appends the following to the data base: PI, QINF, MACH, TINT, 

VAPl, VAP2, RH, RHO, VISC, VINF, REYN. 

RTAT Force Data 

RTAT expects the Keyword NBAL to be available in the data base. If 
NBAL is not found, or if it is less than 1, the entire NBAL loop will be 
skipped. Since the model attitude calculations are part of the NBAL loop, 
NBAL is often set to 1 even in the absence of a balance. 


Extra Equations Before Force 

RTAT has the capability to execute extra equations before the force 
computations. The input specifications for this capability is given in 
APPENDIX 0 along with the algorithms. 


Corrections for Balance Interactions (Sse reference 3) 


RTAT expects to find the following uncorrected or indicated balance 
components in the data base: NF, AF, PM, RM, YM, SF. (See Figure G-1.) These 

components are referred to as "delta" components because they are relative to 
a wind-off zero recording of initial loads. These components are treated 
as a 6 X 1 column vector denoted by [fu], where: 


[fu] - 


AF 

SF 

NF 

RM 

PM 

YM 
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Let the con5>onent8 corrected for Interactions be denoted [p], vdiere: 


[f] 


AFC 

SFC 

NFC 

RMC 

PMC 

YMC 


First Order Interactions 

For first order interactions, the balance calibration establishes the 
following matrix reiationship bet%men correct and indicated '^ponents: 

[fu] - [ci][f] 

where [ci] is a 6 x 6 matrix which is the normalised first order interaction 
coefficient matrix with main diagonal elements of unity. 

Provided [ci] is non-singular, the correct delta components are 
found from: 

[f] - [cii 3 [fu] 

where [cil] is the inverse of [ci]. 
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Second Order Interactions 

For second orier Interactions, the balance calibration establishes the 
following Rtatrix relationship between correct and indicated cong>onentst 


[fu] » [ci][f] + [C2][f2] 


where [C2] is a 6 x 21 matrix which is the normalized second order 
interaction coefficient matrix and [f 2] is a 21 x 1 matrix of product 
combinations of [f] as follows: 


[F2l 


AFC*AFC 

JVFC*SFC 

AFC*NFC 

AFC*RMC 

AFC*PMC 

AFC*YMC 

SFC*SPC 

SFC^NFC 

SFC*RMC 

SFC»PMC 

SFC*YMC 

NFC*NFC 

NFC»RHC 

NPC*PMC 

NFC-YMC 

PMC*BHC 

RNC*PHC 

mC*YNC 

PMC*PMC 

PMC*YHC 

YMC*YMC 
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Provided [ci] is non-singular, the correct delta components eire found 
from: 


[f] = [c1i][fu] - [C1IC2][f2] 

vrtiere [C1IC2] is the product of the inverse of [ci] and [C2], that 
is [ciIC2] = [c1i3[c 2]. This equation must be solved iteratively because 
[f] is expressed in t«»ms of [f 2] which is itself a function of [f]. 

RTAT expects to find the matrices [cil] and [C1IC2] in the data base 
following the Keyword INTR as described in Appendix C. 

Translation for Initial Loads 

Since second order interactions are nonline^u:, the tunnel force emd 
moment ccxnponents must be related to the same origin as the bal^mce calibration. 
Typically, the balance calibration establishes an origin with zero output 
corresponding to zero load. Typically, the tunnel establishes cm origin 
corresponding to an initial load where the initial load is defined as the 
correc'. balance con^onents due to model weight computed for a wind-off zero 
recording. Since the tvmnel subtracts the initial loads from subsequent 
data, the tunnel establishes an origin of zero output corresponding to the 
initial loads. Therefore, a tremslation of axes for initial loads may 
be necessary before the tunnel can use the matrix relationship established 
by the balance calibration. 



APPENDIX G 


Let [fu^] 

[f ] 

o 

[fut] 

[ft] 


denote the indicated initial loads, 
denote the correct initial loads, 
denote the indicated total loads, and 
denote the correct total loads, 


where 


[ft] 


AFTO 

SFTO 

NFTO 

RMTO 

PMTO 

YMTO 


By definition 


[fut] = [pu] + [f ] 

O 


and 


[ft] = [p] + [p^] 


Note that 




WAF 

sin a 

o 







WSF 

din ^ 

o 

cos 

a 

o 





-WNF 

cos (b 
' 0 

cos 

a 

o 





WZRM 

sin (b 
^o 

cos 

a 

o 


WYRM 

cos 


WZPM 

sin a 

o 

- WXPM 

cos d) 
^o 

cos 

a 

o 

WXYM 

sin d) 
^o 

cos 

a 

o 

+ 

WYYM 

sin 

a 

o 
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The second order interaction relationship between correct and indicated 
components for the tunnel is then: 


[ft] = [c1i][fut] - [c1IC2][f2t] 

where [f2t] is just [f 2] based on [ft]. This equation must be solved 

iteratively. 

Note that for initial loads, this becomes 

[f ] = [cil][PU ] - [c1IC2][f 2 ] 
o o o 

where [f 2 ] is just [f 2] based on [f ]. 

o o 

Now, for the i'^ iteration, let 

[e^] = [ciic2][f2t] 

and specifically, for initial loads, let 
[e^] » [C1IC2][F2J 
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The iteration technique is then given by: 


Iteration 

[ft] Approximation 

Error 

0 

Ccii][Fu] + CciiDCfUq^ 


1 

[cii][fu] + [cii][fu^] - [ej 


2 

[cii][fu] + [C1I][PU^] - [Cj^] 

[cp-Ce,] 

3 

[cii][Fu] + [cii][fu ] - [e_] 

O ^ 



This iteration is continued until is less them the specified 

accuracy for all components which is obtained from the balance calibration. 

The correct delta components are then calculated from 

[f] = [ft] - [f^] 

For a wind-off zero, RTAT saves [f^] in a COMMON area. 

RTAT places the following in the data base: AFC, SFC, NFC, RMC, PMC, 

YMC, AFTO, SFTO, NFTO, RMTO, PMTO, YMTO. 

Computation of Sting Deflections 

Sting deflection (or bending) occurs due to loads applied through 
the model. Devices used to measure angles defining balance amd model attitude 
may be so located that they do not record these deflections. Consequently, 
sting deflection angles must be confuted as a function of correct balance 


251 





r * . . 





APPENDIX G 


i. 

I 

1 


loads. The deflection angles are assumed to be small enough so that the 
sting responds elastically, allowing the angle to be described as a spring 
constant times a load. 

RTAT expects the Keyword KDFL to oe available in the data base. If 
KDFL is not found in the data base, it is set to TOTAL. If KDFL is equal to 
DELTA, RTAT computes deflections based on currect delta loads as; 

YAWS = SFC * SFDF + YMC * YI4DF 
ALPS = NFC * NFDF + PMC * PMDF 
PHIS = RMC * RMDF 

If KDFL is equal to DELTA, RTAT next computes the initial deflections as: 
YAWI =0.0 
ALPI =0.0 
PHII = 0.0 

If KDFL is equal to TOTAL, RTAT computes deflections based on correct total 
loads as: 

YAWS = SFTO * SFDF + YMTO * YMDF 

ALPS = NFTC * NFDF + PMTO * PMDF 

PHIS = RMTO * RMDF 

If KDFL is equal to TOTAL, RTAT next computes the negative of the initial 
deflection as: 

YAWI = -YAWS 

ALPI = -ALPS 

PHII = -PHIS 



\ 

j 

! 
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Computation of Balance Attitude 

The attitude of the balance with respect to gravity is determined on the 
basis of a specified input rotation scheme which consists of an ordered 
set of orthogonal Eulerian transformations. Each transformation rotates 
the components of a vector through a specified angle about a specified 
axis. The final result is the transformation of the con 5 )onents of a vector 
from the gravity axis system to the balance axis system. The axis systems 
used are right-hand Cartesian systems. 

Consider a rotation angle denoted by y. The orthogonal trcuisformation 
matrices describing a rotation through amgle Y about each of the three 
possible axes are: 

r 

10 0 
0 cos Y -sin Y 

0 sin Y cos Y 


Roll: 


Rx(Y) = 


Pitch: [Ry(Y)] " 


cos Y 0 
0 1 
sin Y 0 


-sin Y 
0 

cos Y 


Yaw: [Rj,(Y)] = 


cos Y -sin Y 0 
sin Y cos Y 0 
0 0 1 
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The resultant transformation matrix [r] which describes the attitude 
of the balance with respect to gravity then follows from successive applications 
of the appropriate individual transformations. Assuming n rotations, this 
can be written 



] = )][«<>' 

an a 


n-1 


)] . . . t\(Y2)][R3(Y^)] 


where a is either x, y, or z as specified for each y. 

The balance attitude can also be summarized as a single yaw 

rotation YAWG, followed by a single pitch rotation ALPG, followed by a single 
roll rotation PHIG. That is: 


[Rgb] = [Rj^{PHIG)][Ry(ALPG)][R^(YAWG)] 

Substituting and carrying out the indicated multiplications yields: 



Rll 

R12 

R13 

R21 

R22 

R23 

R31 

R32 

R33 


cosa cosU) 

g g 


cos4 sini|< -sin(j) sina cost!) 

g g g g g 


sindi sin<l^ +cos(b sina cos\b 

g g g g ^g 


cosa sinti) 

g g 


cos(() cosijj +sindi sina sinlli 

g g g g ^g 


sincj) cosip -cosd) sina sint^ 

g g g g g 


- sina 


- sin<j) cosa 

g g 


cosdi cosa 

g g 
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Solving for YAWG, ALPG, auid PHIG gives: 


YAWG = arctan (-R12/R11) + ip = -arctan (R12/R11) + ip 

<J q 


ALPG = arcsin (-R13) = -arcsin(R13) 


PHIG = arctan (-R23/R33) + = -arct6m (R23/R33) + (|)^ 


under the following restrictions : 


- 2 - ^^-2 


If 

Rll 

> 

0 

and 

R12 

> 

0, 

*q- 

0 

If 

Rll 

> 

0 

and 

R12 

< 

0, 

♦q- 

0 

If 

Rll 

> 

0 

and 

R12 

= 

0, 

YAWG 

= 0 

If 

Rll 

< 

0 

and 

R12 

< 

0, 


-180 

If 

Rll 

< 

0 

and 

R12 

< 

0, 


180 

If 

Rll 

< 

0 

and 

R12 


0, 


±180 

If 

Rll 

= 

0 

and 

R12 

> 

0, 

YAWG 

= -90 

If 

Rll 

= 

0 

emd 

R12 

< 

0, 

YAWG 

= 90 

If 

Rll 

= 

0 

and 

R12 

= 

0, 

YAWG 

* 0 

If 

Rll 

> 

0 

and 

RI2 

> 

0, 


0 

If 

R33 

> 

0 

and 

R23 

< 

0, 

q 

0 

If 

R33 

> 

0 

and 

R23 

s 

0, 

PHIG 

* 0 

If 

R33 

< 

0 

emd 

R23 


0, 

♦q- 

-180 

If 

R33 

< 

0 

and 

R23 

ss 

0, 

*q- 

180 

If 

R33 

< 

0 

and 

R23 

S5 

0, 

PHIG 

= ±180 

If 

R33 

S 

0 

and 

R23 

> 

0, 

PHIG 

» -90 

If 

R33 

- 

0 

and 

R23 

< 

0, 

PHIG 

. 90 

If 

R33 

m 

0 

and 

R23 


0, 

PHIG 

- 0 
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For the gravity to balance matrix RTAT expects the Keyword NGB 

to be available in the data base followed by the rotation scheme. If NGB 
is not found, RTAT assume a single pitch rotation of zero degrees. The 
elements of the balance attitude transformation cure computed from NGB 
successive applications of the individual transforms. 

RTAT appends the following to the data base: Rll, R21, R31, R12, 

R22, R32, R13, R23, R33, YAWG, ALPG, and PHIG. For the special case of the 
balance attitude with respect to gravity during a wind-off zero recording, 
the angles YAWG, ALPG, and PHIG are given the special neunes YAHZ, ALPZ, 
cuid PHIZ and appended to the data base. 

Height Tare Computations 

RTAT iteratively ccxnputes the nine weight tare factors (three weights 
and six weights times eunns) used in the reduction of force balance data. 

These factors are computed based on wind-off weight tare recordings made at 
various balance attitudes. 

Note that, for the wind-off zero recording, weight tares are referred 
to as initial loads emd are subtracted from each wind-on recording vAiich 
removes the effects of initial loads. Since data 2 ure recorded at attitudes 
other than the wind-off zero attitude, the data must be corrected for delta 
weight tares. 

The model weight is a vector H. Consider the gravity axis system, it is 
apparent that there are no horizontal components of weight and that the 
vertical component is di.rected downward. The weight vector in the gravity 
axis system W is then 
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= (W W , W ) = (0,0, - W) 
g gx, gy' gz 


where W is the magnitude of at attitude [r]. Because the x and y 

components of are zero, only the third column of needs to be 

considered. 


Force Weight Tare Ccxnputations 

Consider tiie third column of [r ^], denoted [rs], vdiich is 

gb 


[r3] = 


- sin ALPG 


- sin a 

g 

- sin PHIG cos ALPG 

= 

- sin (p cos a 

g g 

cos PHIG cos ALPG 


cos cos a 

_ 9 9 _ 


For the wind-off zero recording [r 3] becomes s 




- sin ALPZ 


- sin a 

o 

- sin PHIZ cos ALPZ 

= 

- sin (b cos a 
o o 

cos PHIZ cos ALPZ 


cos cos a 

_ o o_ 


The delta weight tares are defined as the change in balance components 
due solely to model weight W, relative to initial loads, and are computed i 


— — 


— — ■ 

AFTA 


0 

SFTA 

at 

0 

NFTA 


-w 

— — 






r T 

0 

0 

-w 
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where [r ] is [r . ] for the wind-off zero recording. Now let 
o go 



VI 


(sin a - sin a ) 
g o 

[v] = 

V2 

= [R3] - [r3^] = 

(sin (|ig cos ttg - sin <p^ cos a^) 


V3 


- (cos d cos a - cos di cos a ) 


— 


g g ’^o o 


Then the delta force weight tares can be e^ressed as: 


AFTA 


W(sin a - sin a ) 
g o 

SFTA 

» w[v] » 

w(sin ()) cos a - sin d cos a ) 
9 9 o o 

NFTA 


- w(cos cos a - cos d cos a ) 

L - 


g g o o 


Now, recognizing that the force beams of a balance may each sense a 
different portion of the balance's, own weight as well as the model weight, 
this equation may be rewritten as: 


AFTA 


WAF (sin a - sin a ) 
g o 


WAF*V1 

SFTA 

9 

WSF(sin ^ cos a - sin ^ cos a ) 
g g ’^o o 

- 

WSF*V2 

NFTA 


- viNF (cos cos a - cos cos a ) 

_ 9 9 0 o_ 


WNF*V3 
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The solution of this equation may be written 


WAF = 


WSF = 


WNF = 


provided the denominators are not zero. 

In actual practice, several weight tare recordings are made at 
different balance attitudes which results in an overly defined system of 
equations. The solution for the force weight tares will be written in the 
following particular form for consistency with the moment weight tares. 
Assuming n weight tare recordings, the solution is written 


E (AFTA *V1.) 
i-1 ^ 


E (SFTA.*V2.) 
i-1 


E (NPTA. *V3 . ) 
i-1 ^ 
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where the denominators are: 


p 

DTAF = E (VI *V1 ) 

4—1 ^ * 


n 

DTSF = E (V2 *V2 ) 
i=l ^ ^ 


n 

DTNF « E (V3 . *V3 . ) 
i=l " " 

Let DELW represent the delta weight from the balance calibration and let 
D2TF represent a tolerance level for the force weight tares. 

Con«)utationally WAF, WSF and WNF are initialized to zero and the fore, 
weight tares calculated as follows: 


If DTAF > DETF, 


If DTSF > DETF, 


WAF - 


WSF - 


n 


E 

i-1 

(AFTA^*V1^) 


DTAF 

n 


E 

i-1 

(SFTA^*V2^) 


DTSF 

n 


Z 

(NFTA^*V3^) 


If DTNF > DETF, WNF 


DTNP 
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If DTAF < DETF and j WSF ji 

^WNF ^ 


0, WAF = WSF 
0, WAF ■ WNF 


DELW 

DELW 


If DTSF < DETF and 


A»AF 3i< 0, WSF 
^WNF y 0, WSF 


WAF + DELW 
WNF 


If DTNF < DETF and 



¥ 0 , 

¥ 0 , 


WNF = WAF + DELW 
WNF = WSF 


Moment Weight Tare Computations 

If transfer distances X, Y, Z are measured in the balance axis system 
from the bal 2 uice moment center to the model center of gravity* positive in the 
directions of positive thrust* side force, and normal force respectively* 
then the delta moment weight tarn's are obtained by transferring moments as follows: 


1 — ^ 


— 

RMTA 


SFTA*Z - NFTA*Y 

PMTA 

s 

AFTA*Z + NFTA*X 

YMTA 


SFTA*X + AFTA*Y 


SubsLxtuting for AFTA* SFTA* and NFTA gives: 




~ _ 

FMTA 


WSF*V2*Z - WNF*V3*Y 

PMTA 

9 

WAP*V1*Z + WNF*V3*X 

YCTA 


WSF*V2*X + WAF*V1*Y 

_ _ 

1 

L_ J 
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which can be rewritten as 


RWTA 

PMTA 

YMTA 


WZRM*V2 - WYRM*V3 
WZPM*V1 + WXPM*V3 
WXYM^V2 + WYYM*V1 


r 


Writing the individual equations gives: 


[rMTa] = [V2 - V3] 


WZRM 

WYRM 


[pmTA] = [Vl V3] 


WZPM 

WXPM 


[yMTA] = [V2 Vl] 


WXYM 

WYYM 


Each of these equations is solved for the weight tare factors by first premultiplying 
both sides of the equation by the transpose of the row vector on the right 
containing the V terms. The solution then involves inverting the matrix 
obtained from the transpose of the row vector times the row vector itself. 
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WYYM 


rVMTA *\T0 \ 

Z (V2.*V2 )* Z (YMTA.m.) » Z (V2.*Vl.)* Z 
i=l ^ i=l ^ ^ i=l ^ ^ i*l 


DTYM 


where the denominators are: 


n n n n 

DTRM = Z (V2.*V2.)* 2 (V3.^V3J - Z (V2.*V3J* Z (V2,*V3.) 

.. 1 X.,1 1 .,1 l.'l 1 

1=1 1=1 1=1 i=l 


n n 

DTPM = I (V1.*V1,)* I (V3.*V3.) 

. - 1 i . n ^ 1 

1=1 1=1 


n n 

Z (V1.*V3.)* Z (V1.*V3 ) 
i=l ^ ^ i=l ^ ^ 


n n n n 

DTYM = I (V1.*V1.)* E (V2.*V2.) - E (V1.*V2.)* E (Vl . *V2 ) 
i=l ^ ^ i=l ^ ^ i=l ^ ^ i=l ^ ^ 


Let DETM represent a tolerance level for the moment weight tares. 
Computationally, WZRM, WYRM, WZPM, WXPM, WXYM, and WYYM are calculated 
as follows: 
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If IdTRmI > DETM, 


n n 

^ (V3.*V3.)* Z (rmTA.*V2.' 
1-1 ^ ^ • -1 i 1 

WZRM = — ^ — 

DTRM 


n 


n 


^ (V2.*V3.)* ^ (RMTA.*V3.) 
1=^1 " ^ i=l " ^ 


and 


WYRM = 


n n 

S (V2.*V3.)* Z (RMTA.*V2.) 
i=l ^ ^ i=l ^ ^ 


DTRM 


n n 

(V2 *V2.)* I (RMTA.*V3.) 

_-i ^ 1 11 


If |dTPm| > DETM, 


n n 

S (V3 *V3.)* L (PMTA.*V1.) 
i =1 ^ ^ ' —1 ^ ^ 

WZPM = — 

DTPM 


n n 

I (VI *V3J* Z (PMTA.*V3.) 
-1 ^ 1 11 


and 


WXPM 


n n 

Z (V1.*V1.)* Z 
-ill . 


(PMTA^*V3^) 


DTPM 


n n 

E (VI *V3.)* I VI ,"A.*V1.) 
i=l ^ i*l ^ ^ 
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If IdTYmI < DETM 


WXYM 


n n n 

^ (V1.*V1.)* Z (YMTA.*V2.) - Z (V2.m.)* Z (YMTA.*V1 ) 

i=l ^ ^ i=l ^ ^ i=l ^ ^ i=l ^ ^ 


DTYM 


and 


WYYM 


n n n n 

S (V2.*V2.)* Z (YMTA.*V1.) - Z (V2.*V1.)* Z (YMTA.*V2.) 
1-1 " * 1-1 ^ * 1.1 " ‘ i-i ‘ ‘ 


DTYM 


If |dTRm| < DETM, WZRM = 0 and WYRM = 0 


If 1 dTPm| < DETM, WZPM = 0 and WXPM = 0 


If |dTYm| < DETM, WXYM = 0 and WYYM = 0 


If |dTRm 1 < DETM and 


If |dTPm| < DETM and 


If |dTYm| < DETM and 


! |dTPm| > DETM, 
I DTYM I > DETM, 

f iDTYMj > DETM, 
|dTRM| > DETM, 

! |dTKM| > DETM, 
|dTPm| > DETM, 


WZRM = WZPM 
WYRM = WYYM 


WXPM = WXYM 
WZPM s WZRM 

WYYM * WYRM 
WXYM = WXPM 
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The weight tare moment arms can then be calculated (for reference 
purposes only) as follows 


YRM 


WYRM 

WNF 


ZRM 


WZRM 

WSF 


XPM 


WXPM 

WNF 


ZPM 


WZPM 

WAF 


XYM 


WXYM 

WSF 


YYM = 


WYYM 

WAF 
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4 


# 

r 


RTAT computes the weight tare factors based upon correct balance 
components, which in turn are based upon assumed initial loads equal to the 
wind-off zero loads. The newly computed weight tare factors are thsn used to 
recompute the initial loads, which are then compared to the original initial 
loads. If the new initial loads are sufficiently close to the old initial 
loads, they are assumed to have converged, and the tare computations are 
accepted. If the new initial loads are significantly different, the old 
initial loads are replaced with the new initial loads, which are then used to 
recompute correct balance components to start another iteration. This 
process continues until the initial loads converge or for a maximum of five 
iterations. The accuracy required for each component is obtained from the 
balance calibration. 

HTAT saves in a COMMON area the nine weight tare factors, the six 
initial loads, and the wind-off zero attitude. 


Correction for Weight Tares 

RTAT obtains the nine weight tare factors, the six initial loads, and the 

wind-off zero attitude from a COMMON area. 

For the balance at any attitude [R . ]» [pTARE] is computed from the 

gb 

previously given equations v^ere: 


[ftare] 


AFTA 

SFTA 

NFTA 

RMTA 

PMTA 

YMTA 


268 









APPENDIX G 


Let [fBAL] 


denote the earodynamic loads in the balance axis system 


where 


[fbal] 


AFBA 

SFBA 

NFBA 

RMBA 

PMBA 

YMBA 

— — 


The aerodynamic loads in the balance axis system are computed by 
subtracting the weight tares from the correct delta balance loads: 

[fbal] = [f] - [ftare] 

RTAT appends the following to the data base: VI, V2, V3, AFTA, 

SFTA, NFTA, RMTA, PMTA, YMTA, AFBA, SFBA, NFBA, RMBA, PMBA, YMBA. 

Computation of Model Attitude 

The gravity axis to balamce axis transformation been 

fully developed. Similar transformation matrices from wind axis to gravity 

axis [r ] and from balance axis to model axis [r. ] may be defined on the 
wg Dm 

basis of specified input rotation schemes. 

The wind to gravity matrix [r ] describes the tvmnel flow angularity 

wg 

and is defined as an upflow or pitch rotation of magnitude ALPU and a crossflow 

or yaw rotation of magnitude XFLO. As presently constituted [r ] does not 

wg 

include a roll rotation term, [r ] is built into RTAT. 

wg 
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For the balcince to model matrix RTAT aspects to find tlie Keyword 

NBM in the data base followed by the rotation scheme. If NBN is not founds 
RTAT assumes a single pitch rotation of zero degrees. The elements of the 
balance to model transformation are computed frcm NBM successive applications 
of the individual transforms. 

The attitude of the model with respect to the wind can then be described 
by the transformation matrix [^^3 given by: 

The model attitude also be summeurized as a single roll 

rotation PHIW, followed by a single yaw rotation YAWW, followed by a single 
pitch rotation ALPW. That is: 


[R^^] = [ Ry (ALPW) ][r^(YAWW)][r^ (PHIW)] = [Ry (a^) ] [R,, ( V ^^^x ^ 


Substituting the elementary transforms and carrying out the indicated 
multiplications yields: 



WMll 

WM12 

WM13 

WM21 

WM22 

WM23 

WM31 

WM32 

WM33 
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cosa cosii^ 

w w 


-cosa sinij^ cosi|) -sina sin(J> cosa sinit sin(b -sina cosd> 

WWW W ^W WWW w w 


sinip 


w 


sina cosiL 
w w 


COSl(^^COS(()^ 


-sina sinit cosd) +cosa sind) 
WWW w w 


-cosiJ^^sin(|>^ 


sina sinil^ sind) +cosa cosd) 
w ^w w w ^w 


Solving for YAWW^ ALPW, and PHIW gives: 


PHIW = -arctan(WM23/WM22) + <j> 


YAWW = arcsin (WM21) 


ALPW = arctan (WM31/WM11) + a 


under the following restrictions: 


7T 7T 

- y £ YAWW < y 
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If WM22 > 0 and WM23 >0, * 0 

If WM22 > 0 and WM23 <0, = 0 

If WM22 > 0 and WM23 = 0, PHIW = 0 

If WM22 < 0 and WM23 >0, <1)^ = -180 

If WM22 < 0 and WM23 <0, 0^ = 180 

If WM22 < 0 and WM23 = 0, PHIW = ±180 

If WM22 - 0 and WM23 > 0, PHIW = -90 

If WM22 = 0 and WM23 < 0, PHIW = 90 

If WM22 = 0 and WM23 ® 0, PHIW = 0 

If WMll > 0 and WM31 >0, = 0 

If WMll > 0 and WM31 <0, * 0 

If WMll > 0 and WM31 » 0, ALPW « 0 

If WMll < 0 and WM31 >0, * -180 

If WMll < 0 and WM31 <0, - 180 

If WMll < 0 and WM31 » 0, ALPW = ±180 

If WMll = 0 and WM31 > 0, ALPW = -90 

If WMll * 0 and WM31 < 0, ALPW * 90 

If WMll - 0 and WM31 » 0, ALPW » 0 
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The anale of sideslip is computed from the definition: 


BETA » -YAWW 

RTAT appends the follo«/ing to the data base: WMlli WM21, WM31, WM12, 

WM33, WM32, WM13, WM23, WM33, PHIH, YAWW, ALPW, BETA. 

Computation of Components in Model Axis System 
RTAT obtains the correct balance axis components from the data base 
and makes use of the rotation matrix [r ^ ] to rotate them to the model axis 
system: 


"*AFMA^“ 


AFBA ’ 

SFMA^ 


SFBA 

NFMA^ 


NFBA 

^ 


^ , 




-RMBA 

PMMA^ 

■ [^3 

PMBA 

-YMMA^ 


-YMBA 



1^ ^ 
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The minus signs in the moment equation resiU-t from the inherent opposite- 
handedness of the conventional wind tunnel balance force and moment sign 
conventions. The subscripts on the model axis con^nents are for clarity in 
the following presentation. 

RTAT obtains the transfer distances XBAR, YBAR, and ZBAR from the data 
base. These describe the transfer from the balance moment center to the model 
moment reference center. The tremsfcr equations are given by: 


AFMA^ 


"aFMAj^ 

SFMA^ 


SFMAj^ 

NFMA^ 


NFMAj^ 



RMMAj^ + NFMAj^*YBAR - SFMAj^*ZBAR 

PMMA^ 


PMMAj^ - NFMAj^*XBAR - AFMAj^*ZBAR 

YMMA^ 


YMMAj^ - SFMAj*XBAR - AFMAj^*YBAR 

— -J 


_ — 


The subscripts on the model eucis components are for clarity in the 
following presentation. 

RTAT places AFMA, SFNA, NFMA, RMMA, PMMA, and YMMA in the data 

base. 
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Blockage and Jet Boundary Corrections 
The blockage corrections are based on reference 4 and the jet boundary 
corrections are based on reference 5. 

RTAT expects to find the Keyword BLK in the data base . BLK is a flag 
used to control the computation of both blockage and jet boundary corrections. 

A value Oi. zero omits the corrections while a value of unity applies the 
corrections . 

RTAT expects to find the following variables in the data base: KWI, 

KBI, J2, J3, PI, QINF, MACH, REYN, RHO, VINF, B, S, ALPW, BETA, AFMA, 

SFMA, I “MA. 

For convenience, define the tunnel :-araweter8 uncorrected for blockage as 


PIPR « PI 
QPR - QINF 
MPR « MACH 
RHT , 

‘r 1 

VPR - VINF 

The apparent model axis forces are then ccs^uted: 


ACAN • AFNA/<^PR*S) 
ACYN - 8FMA/(QPR*S) 
ACNM ■ NFHA/(QPR*S) 
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These apparent force coefficients are then rotated to the stability and 
wind axes: 


’acds 


"acam 

ACYS 

= [r^(-ALPW)] 

ACYM 

ACLS 


ACNM 


ACDW 


ACDS 

ACYW 

= [Rj^(BETA)] 

ACYS 

ACLW 

L 

ACLS 


The corrected apparent drag coefficient is then calculated; 


CDPR = ACDW - 


ACLW*ACLW^S 


/ B 

n (jj) 


The blockage corr :tion factor is then 


7^1 t KBI 
(1 - XM2)^/^ 


1 + 0.4*XM2 

+ ( ^)*CDPR*KJ 

1 - XM2 


where 


XM2 B MPR*MPR 


and 


KJ 


4*C 


1 

7x10 
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where the cross sectional area of the 7- by 10-foot high speed 

tunnel test section. (See reference 1.) 

The tunnel parameters are then corrected for blockage: 

PI =• PIPR*(1 - 1.4*XM2*K) 

QINF = QPR*(1 + (2 - XM2)*K) 

MACH = MPR*(1 + (L + 0.2*XM2)*K) 

RHO = R0PR*(1 - XM2*K) 

VINF » VPR*(1 + K) 

REYN = RNPR*(1 +(1 - 0.7*XM2)*K) 

Hie jet boundary corrections require the calculation of the apparent 
lift coefficient: 

j : = ACLW*{2~) 

The jet boundary correction factors for angle of attack and pitching 
moment are then: 
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DELA = CLC*J2 


DELM = CLC*J3*QINF*S*CBAR 


The angle of atte >. eind model eucis pitching moment are then corrected 
for jet boxindary effect: 


ALPW = ALPW + DELA 
PMMA^ = PMMA2 + DELM 

The subscripts on PMMA are for clarity in the following presentation. 

RTAT updates the data base values of the following varied)les: PI, 

QINF, MACH, RHO, VINE, REYN, ALPW, PMMA. RTAT appends the following variables 

to the data base; PIPR, QPR, MPR, RNPR, ROPR, VPR, CDPR, CDC, DELA, DELM. 

Base and Chamber Pressure Corrections 

For base pressure, RTAT expects the Keyword NBAS to be available in the 
oata base. If NBAS is not found or if it is less than 1, no base pressure 
computations will be performed. RTAT expects the base pressure specifications 
to follow NBAS in the data base. RTAT obtains from the data base the vciriables 
S, PI, QINF, and ALPW. 

Letting PB^ denote the i^^ base pressure, the individual pressure 
coefficients cure computed as: 

PB, - PI 

CPBi » — = 

QINF 
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The forces and moments to be applied as a correction are computed as: 


AFB 

SFB 

MFB 

RMB 

?21B 

YMB 


Z(PB. “ PI) 
i 


-Area 


AF. 


Area 


SF. 


Area^^ 

NFi 

(Area*Arm) 


(Area*Arm) 


RM. , 
1 


PM. 


(Area*Arm) 




where the summation i is over all base pressures for which the correction flag 
is turned on and where the column vector on the right side and the correction 
flag are obtained from thr input specifications. 

The forces and moments to be computed but not to be applied as a 
correction are given by: 


xafb" 


-Area.- 

AF. 

1 

XSFB 


Area__ 

SFi 

xnfb 

= Z(PB. - PI) 

1 

XRMB 

1 

1 

(Area*Ann) 

RMi 

XPMB 


(Area*Arm)_„ 

FM , 
X 

XYMB 


(Area*Arm) ^ 
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where the summation i is over all base pressures for which the correction 
flag is turned off and where the column vector on the right side and the 
correction flag are obtained from the input specification. 

The axial force terms are also expressed in terms of axial and drag 
coefficients as: 

CAB = AFB/(QINF*S) 

CDB = CAB*COS (ALPW) 

XCAB = XAFB/(QINF*S) 

XCDB = XCAB*COS (ALPW) 

RTAT appends the following to the data base; CPBi, AFB, SFB, NFB, RMB, 

PMB, YMB, XAFB, XSFB, XNFB, XRMB, XPMB, XYMB, CAB, CDB, XCAB, XCDB. 

For chamber pressures, RTAT expects the Keyword NCBR to be available 
in the data base. If NCBR is not found, or if it is less than 1, no chamber 
pressure ccanputations will be performed. RTAT expects the chamber pressure 
specifications to follow NCBR in the data base. RTAT obtains from the data 

base the variables S, PI, QINF, and ALPW, 

Letting PC^ denote the i^^ ch 2 unber pressure, the individual presstire 

coefficients are computed as; 

PC. - PI 

CPCi = — 

QINF 
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The forces and moments to be applied as a correction are computed as: 


4 


liFCH 


-Area,„ 

AF. 

1 

SFCH 


Area„_ 

SF. 

1 

HFCH 

= j:(pc. - PI) 

Area.— 

NF. 

1 

RMCH 

X 

1 

(Area^Arm) 

PMCH 


(Area*Arm)_„ 
PM , 
X 

YMCH 


(Area*Am)^ 
YM . 

_ X_ 


where the summation i is over all chamber pressures fcr which the correction 
flag is turned on and where the column vector on the right side and the correction 
flag are obtained from the input specifications. 

The forces and moments to be computed but not to be applied as a correction 
are given by: 


XAFC 


-Area,„ 

AF. 

1 

XSFC 


Area_„ 

SF. 

1 

XNFC 

= I(PC^ - PI) 

i 

Area.„ 

NF. 

1 

XRMC 

(Area*Arm)jy^ 

*'"i 

XPMC 


(Area*Arm)_„ 
PM , 
1 

XYMC 


(Area*Arm) 

YMl 


where the summation i is over all chaunber pressures for which the correction 
flag is turned off and where the column vector on the right side and the 
correction flag are obtained from the input specifications, 
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The axial force terjns are also expressed in terms of axial and drag 
coefficients as: 

CAC = AFC/(QINF*S) 

CDC = CAC*COS (ALPW) 

XCAC = XAFC/(QINF*S) 

XCDC = CDC^COS (ALPW) 

RTAT appends the following to the data base: CPCi, AFCH, SFCH, NFCH, 

RMCH, PMCH, YMCH, XAFC, XSFC, XNFC, XKMC, XPMC, XYMC, CAC, CDC, XCAC, XCDC. 

Computation of Model, Stability, and Wind Axis Components 

RTAT obtains the following from the data base: AFMA, SFMA, NFMA, 

RMMA, PMMA, YMMA, AFB, SFB, NFB, RMB, PMB, YMB, AFCH, SFCH, NFCH, RMCH, 

PMCH, YMCH, ALPW, BETA, QINF, S, B, CBAR. 

The forces and moments are corrected for base and chamber pressures 
to give the final corrected model axis components: 


AFMA~ 


^MA^ 




^CH 

SFMA 


SFMA^ 


SFB 


SFCH 

NFMA 


NFMA 2 


NFB 


NFCH 


= 


- 


- 


RMMA 


RMMA^ 


RMB 


RMCH 

PMMA 


PMMA^ 


PMB 


PMCH 

YMMA 


YMMA^ 

1 

YMB 


YMCH 


i 


1 j 
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.'he stability axis components are obtained by rotating the model cucis 
components through minus the angle of attack: 


DRAG 



SFSA 

= [r (-ALPW)] 

SFMA 


y 


LIFT 


NFMA 


~rmsa' 


-RMMA ’ 

PMSA 

= [Ry(-ALPW)] 

PMMA 

-YMSA 


-YMMA 


The wind axis components are obtained by rotating the stability eucis 
components through the angle of sideslip: 


AFWA 


DRAG 

SFWA 

= [r (BETA)] 
z 

SFSA 

OTWA 


LIFT 


“^RMti/A 


-RMSA 

PMWA 

= [r (BETA)] 
z 

PMSA 

-YMWA 


-YMSA 
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Computation of 
Let [fma], [fsa], ind [fwa] denote 


Coefficients 

the model, stability, and wind 


axis components: 


[fma] 


AFMA 

SFMA 

NFMA 

RMMA 

PMMA 

YMMA 


[fsa] = 


DRAG 

SFSA 

LIFT 

RMSA 

PMSA 

YMSA 
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[fwa] 


AFWA 

SFWA 

NFWA 

RMWA 

PMWA 

YMWA_ 


Define the 6x6 main diagonal matrix [c] as ; 


1 

QINF*S 

1 

QINF*S 


[c] 


1 

QINF*S 

1 

QINF*S*B 

1 

QINF*S*CBAR 

1 

QINF*S*B 

_J 


where all of the elements off the main diagonal are zero. 

Let [csa], and [cwa] denote the model, sted>ility, and wind 

axis coefficients: 
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The model, steibility, and wind axis coefficients are then computed as: 

[cma] = [c][fma] 

[csa] = [c][fsa] 

[cwa] = [c][fwa] 

The lift to drag ratio and lift squared are also computed as; 

L/D = CL/CD 
CLSQ = CL*CL 

RTAT updates the data base vlaues of the following variables: 

AFHA, SFMA, NFMA, RMMA, PMMA, YMMA. m‘AT appends the following variables 
to the data base: DRAG, SFSA, LIF-T, RMSA, PMSA, YMSA, AFWA, SFWA, NFWA, 

RiWA, PMWA, YMWA, CA, CY, CN, CRM, CM, CYM, CD, CYS, CL, CRMS, CMS, CYMS, 
CDW, CYW, CLW, CRMW, CMW, CYMW, L/D, CLSQ. 
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Pressure Coefficient Arrays 

RTAT expects the Keyword NCP to be available in the oa ^ base. If 

it is not found or if it is less than 1, no pressure coefficient arrays 

will be computed. RTAT obtains PI and QINF from the data base. 

Let P denote the name of the input pressure and C denote the name 

P 

of tht! output pressure coefficient, then the pressure coefficient is computed 
as: 


P. - PI 

C - 

^i QINF 

where i goes from 1 to NSIZE which is obtained from the input specification. 
Note that NSIZE consecutive data items in the data base starting with P 
%d.ll be replaced by the corresponding and thus P will not be available 

for further processing or output. 

Pressure Ratio Arrays 

RTAT expects the Keywrii NRTO to be available in the data base. If 
it is not found or if it it. less than 1, no pressure ratio arrays will be 
computed. 

Let P denote the nasm of the inp'it pressure and denote the name 
of the output pressure ratio, then the pressure ratio is cce^uted as: 

" Scalar 
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where i goes from 1 to NSIZE which is obtained from the input specification 
as IS the neune Scalar. Note that NSIZE data items will be appended to the 
data base. 


Flowmeter Confutations 

RTAT expects the Keyword NFLO to be available in the data base. If 
NFLO is not found or is less than 1, r\f ilowmeter computations will be 
performed. RTAT obtains the following variables from the data base: HI, 

TT, FPn, FDPn, and FTn where n is the flowmeter number. 

RTAT contains complete bles of parameters for all of the venturi- 
type flowmeters available for use in the 7- by 10-foot high speed tunnel. 

The flowmeter fluid viscosity is computed as 

FT 

Z.'Un = ((.1211*10" ) + ((4/3)*10"^*FP))‘ 


The flowmeter diameter ratio is given by: 



APPENDIX G 


1 

I- 

f The velocity of approach factor is computed as: 


Fn= 1/ 



The throat static pressure is computed asi 


i’’t’2n- FPn~ FDPn 


The static pressure ratio is computed as: 


0 

SPRn = FP2n/FPn 


The flowmeter expansion factor is con?)uted as: 


Yn = YAAn + YABn*Rn 


The flowmeter temperature is converted to degrees Fahrenheit 


by: 


TFjj = FTn- 459.688 


The flovnneter supercompressibility factor is computed as: 



Nr 
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B2 = .61113723*10“'^-.67725162*10‘^*TFn+. 2912419*10"®*TFn^-. 707898l5*10"^^*TFn^ 

B3 = -.1 5730387*10"^+. 10578106*10"^*TFn-.14254673*10"^^*TFn+.81196439*10°^^*TFn^ 
B4 = -.76772236*10"^^-. 12788883*10"^*TFn+.59147764*10“^^*TFn^-.36507332*10"^^*TFn^ 
PF = FPn - 25. 

2 3 

SFMn = 1 + B2 * PF + B3 * PF + B4 * PF 

The discharge coefficient, weight flow rate, and Reynolds niimber are 
iteratively computed as: 

WPn = DCn*Fn* */MS(SFMn)* A2n* Yn* 158.1948 * /aBS (FP n*FDPn/FTn) 

RNn = (4 * WPn)/(ZMUn * 3.14159 * D2n/12) 

3 

DCn = E COEF. **RNn^ 
i=0 ^ 

Normalized temperature is computed as; 

THTn = TT/ (459.688 + 59) 

Normalized pressure is computed as: 

LAMn = HI/2116.8 

Normalized weigitt flow rate is coniputed as: 


WPN 


WPn * ( 


✓rariT 

LAMn 


) 
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RTAT appends the following items to the data base: AMUn, Din, D2n, 

BETn, Fn, FP2n, SPRn, Yn, TFn, SFMn, ASQn, DCn, WPn, RNn, THTn, LAMn, and 
WPNn where n is the flowmeter number. 

Jet Exhaust Computations 

When required, RTAT may be assembled with a special routine to permit 
real time calculation of and correction for jet exhaust effects. These 
computations are usually model dependent. 

Unique Computations 

When required, RTAT may be assembled with a special routine to perform 
computations unique to a particular model. 

Extra Equations After Force 

RTAT has the capability to execute extra equations after the force 
computations. The input specifications for this capability is given in 
APPENDIX D along with the algorithms. 

Real Time Displays 

RTAT expects the data item NDSP to be available in the data base. 

If NDSP is not found, or if it is less than 1, only ALPW and BETA will be 
displayed on the Automatic Angle Panels. The input specifications are given 
in APPENDIX C. 

When a code number is entered into the thumbwheels beside one of the 
displays, the data value associated with that code number is displayed on 
that display. The displayed values are updated during every execution of RTAT. 
The exact timing depends on the OAP and RTAT input setup constants, but is 
normally approximately once every second. The displays cannot be updated 
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while the RTAT task is executing a recorded data point, resulting in a 
delay of two to thirty seconds depending on the amount of printout and the 
type of plotting produced. 


Line Printer Output 

RTAT produces two different types of printout: point-by-point 

printout and run-by-run summary printout. 

Point-By-Point Print 

RTAT expects the data item NPG to be available in the data base. If 
NPG is not found, or if it is less than 1, no point-by-point line printer 
output will be produced. 

Point-by-point printout is produced for every recorded data point. 

Run-By-Run Summary Print 

RTAT expects the data item NGP to be available in the data base. If 
NGP is not found, or if it is less than 1, no run summary line printer output 
will be produced. 

The actual run summary printout is produced only on request in response 
to the thumbwheel entry of a specific data identification code. 

The data requested on the specifications is saved on a point-by- 
point basis on files on the disk. The print segment will give a warning 
message before overflow of these files occurs. 

Plot Output 

RTAT produces two different types of plots: point-by-point plots 

and summary plots. 
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Point-By-Point Plots 

RTAT expects the data item NPLT to be available in the data base. If 
NPLT is not found, or if it is less them 1. no plots will be produced. 

If the NPLT specifications are used to plot pressure arrays, then RTAT 
expects the data item MPLT to be available in the data base. If MPLT is 
not found, or if it is less than 1, no additional plots will be produced. 
These specifications may only be used to plot additional pressure arrays. 

For pressure arrays, a complete plot is generated for each test point. 
For force data a complete plot is generated for each run. Hard copies are 
automatically generated as appropriate. 

Summary Plots 

RTAT expects the data item SPLT to be available in the data base. 

If SPLT is not found, or if it is less than 1, no summary plots will be 
generated. 


RTAT Output Tape 

RTAT will write an answer tape in standard interface format. This 
tape may be used for subsequent processing. 
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Figure G-1.- Force and mc»nent model euces with 
ixisitive directions shown. 
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RTAT OPERATING PROCEDURES 
Initial OAP Startup 

Before starting RTAT, it is necessary for the current version of the 
7- by 10-foot tunnel operating system to be running on the Sigma 3 computer. 
An appropriately prepared setup deck should then be loaded into the card 
reader. The OAP will start execution and wait for an Enter on the System 
Control Panel. The following sequence of operations must then be performed: 

(a) On the System Control Panel, set the Data Ident thumbwheels to 
00. Select the ON position for Select Options, List, Tape, Request, Free. 
Verify the proper setup of the other control pemels. 

(b) Press the Enter button on the System Control Panel. OAP will 
read the OAP setup cards and perform a consistency check. OAP will report 
cuiy inconsistencies and request that they be resolved before proceeding 
further. 

(c) Select the ON position for Calibrate on the System Control Patnel 
and press the Enter Button. OAP will perform an analog system calibration 
and report the results. 

(d) Select the ON position for Dau Data on the System Control Panel 
and press the Enter button. Press the Ready/Cycling button on the System 
Control Panel. OAP will record a zero data point and trigger the RTAT tasks. 
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Initial RTAT Startup 

The data ident codes currently recognized by RTAT are presented in 
Figure H~l. 

After GAP triggers the RTAT tasks, perform the following steps: 

(a) On the System Control Panel, set the Data Ident thumbwheels 
to 09 and pressure the Enter button. RTAT will read the RTAT setup cards 
and perform a names record execution and report the results. 

(b) Set the Data Ident thumbwheels to 00 and press the Enter button. 
Press the Ready /Cycling button. GAP will record a zero data point and RTAT 
will compute the zero data point. 

RTAT Setup Change 

Whenever it is desired, a new RTAT setup may be performed. The 
recommended method of accomplishing this is to read in cui entire setup deck 
as follows: 

(a) Place the new RTAT setup deck in the card reader. 

(b) Gn the System Control Panel, set the Data Ident thumbwheels 

to 09 and press the Enter button. RTAT will read the RTAT setup cards and 
perform a names record generation execution and report the results. 

If the only changes required are the addition of new D-type constants 
or changes to the values of C-type constants, it is possible to input only 
the new cards and not an entire setup deck. This may be accomplished by 
using a Data Ident thumbwheel setting of 29 in step (b) above « This method 
is not recommended for anyone without a thorough knowledge of the data base 
because it may lead to subtle errors in the interpretation of the new 
setup by RTAT and may invalidate the real time data reduction. 
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Responding to a Change in GAP Setup 
Certain changes in the GAP setup require a response from RTAT in the 
form of a new names record execution without the need for a change in the RTAT 
setup deck. This will occur, for example, if the length of the GAP write- 
scan table is changed by adding or deleting channels or by turning the 
scanivalves on or off. (Note: If the scanivalves are turned on or off, 

RTAT requires that the GAP PVID bits also be turned on or off correspondingly 
in order for RTAT to correctly interpret the status of GAP. That is, with the 
scanivalves turned off but the PVID bits turned on, RTAT will not generate a 
name for the analog channels indicated by the PVID bits.) 

RTAT will not automatically respond to a change on GAP setup, the operator 
must so direct RTAT. This is accomplished by the following step: 

(a) On the System Control Panel, set the Data Ident thumbwheels to 
19 and press the Enter button. RTAT will perform a names record generation 
execution and report the results. 

This procedure may also be used when there is no change in GAP setup in 
order to obtain an additional names record execution report. 

Wind Off Zero Records 
A wind off zero point is recorded as follows: 

(a) On the System Control Panel, select the ON position of Dau Data, 
set the Data Ident thumbwheels to 00 and press the Enter button. 

(b) Press the Ready/Cycling button. GAP will record a zero data 
point and RTAT will compute a wind off zero data point. 
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Weight Tare Records 

A wind off zero must be recorded prior to recording any weight tare 
data. At least two weight tare points must be recorded at attitudes 
sufficiently different from each other and from the wind-off zero attitude 
to allow proper angle resolution during the weight tare factor computations. 
A normal weight tare run may include data recorded at a half-dozen different 
attitudes. 

A weight tare point is recorded as follows; 

(a) On the System Control Panel, select the ON position of Dau Data# 
set the Data Ident thumbwheels to 0.1 and press the Enter button. 

(b) Press the Keady/Cycling button for each weiyht tare point to be 
recorded. OAP will record a data point and RTAT will compute a weight tare 
point and save it on a special disk file for use in weight tare factor 
computations. 


Weight Tare Factor Computations 

After a set of weight tare points have been recorded# the operator 
must direct RTAT to compute the weight tare factors. This is accomplished 
as follows: 

(a) On the System Control Panel# set the Data Ident thumbwheels to 
21 and press the Enter button. RTAT will read the weight tare data it saved 
on the disk# compute the weight tare factors# generate a report# ^md save the 
weight tare factors in a COMMON area for later use in the conqputations. 
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Wind-On Data Records 

A wind-off zero record must be recorded prior to recording any wind-on 
data. If weight tare corrections are to be applied, the weight tare 
factors must be available to RTAT. 

Wind-on data are recorded as follows: 

(a) On the System Control Panel, select the ON position for Dau Data, 
set the Data Ident thumbwheels to 04, and press the Enter button. 

(b) Press the Ready/Cycling button for each wind-on data point to be 
recorded. OAP will record a data point and RTAT will compute the wind on 
data point. 


Calibration Data Records 
Calibration data are recorded as follows: 

(a) On the System Control Panel, select the ON position for Dau 
Data, set the Data Ident thumbwheels to 08, and press the Enter button. 

(b) Press the Ready/Cycling button for each calibration data point 
to be recorded. OAP will record a data point and RTAT will compute a 
calibration point and save it on a special disk file for use in an interactive 
calibration workup session. 

Interactive Calibration Conqputations 
After a set of calibration points have been recorded, the operator must 
direct RTAT to workup the calibration constants. This is accomplished as 
follows: 

(a) On the System Control Panel, select the ON position for Dau Data, 
set the Data Ident thumbwheels to 28, and press the Enter button. RTAT 
will enter an interactive mode and request inputs from the operator on the 
interactive device. The interactive commands are described in APPENDIX E. 
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Data Base Checkpoint Dumps 

RTAT has the capability to generate a report on the complete contents 
of the data base at several points In the calculation sequence. This 
capability exists for any Data Ident having 0 for Its left digit and is 
accomplished by chemglng that left digit to a 1. 

SuIITm^u:y Output Generation 

RTAT has the capability to save certain information on a special 
disk file. The operator must direct RTAT to process this summary file. This 
is accomplished as follows: 

(a) On the System Control Panel, set the Data Ident thundawheels 

to 24 and press the Enter button. RTAT will read the suunary data from the 

special disk file, generate a sunmary report, and generate summary plots. 

Pressurized Zero Data Records 

A wind'off zero point must be recorded before a pressurized zero is 
recorded. For powered model testing, a pressurized wind^^ff zero is accotplished 
as follows after pressurization: 

(a) On the System Control Panel, select the ON position for Dau Data, 

set the Data Ident thumbwheels to 05, and press the Enter button. 

(b) Press the Ready/Cycling button for each pressurized wind-off 
zero to be recorded. OAP will record a data point and RTAT will compute the 
pressurized data point as a data point and not a zero point (which is what is 
desired) and generate a report. 
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Wind-On Powered Data Records 

For powered model testing, powered wind on data are recorded as follows: 

(a) On the System Control Panel, select the ON position for Dau Data, 
set the Data Ident thumbwheels to 06, and press the Enter button. 

(b) Press the Reudy/Cycling button for each powered data pcint to be 
recorded. OAP will record a data point and RTAT will compute a powered 
data point and generate a report. 

Static Powered Data Records 

For powered models, static powered data are recorded as follows: 

(a) On the System Control Panel, select the ON position for Dau 
Data, set the Data Ident thumbwheels to 07, and press the Ente** button. 

(b) Press the Ready/Cycling button for each static powered data 
point to be recorded. OAP will record a data point. 

Note that RTAT does not currently support static powered data. 

Interactive Static Powered Data Computations 

After a set of static powered calibration points have been recorded, 
the operator must direct RTAT to workup the calibration constants. This is 
accomplished as follows: 

(a) On the System Control Panel, select the ON position for Dau 
Data, set the Data Ident thumbwheels to 27, and press the Enter button. 

Note that RTAT does not currently support interactive torkup of static 
powered calibration data. 
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System Calibrate Cooptation 

When OAP performs an analog system calibration, it vrrites the system 
calibrate information on a special disk file. Additional processing by 
RTAT is required before RTAT can use the system calibrate data. This is 
accomplished as follows: 

(a) On the System Control Panel, select the ON position for Dau 
Data, set the Data Ident thumbwheels to 40, and press the Enter button. 

RTAT will read the OAP system calibration data, process it, generate a 
report, and save it in a form vrtiich will enable RTAT to convert from 
uncorrected counts to corrected millivolts. 

Link to Central Computer Con^lex 
RTAT has the capability to save its inputs and results as a SIF 
on a magnetic tape. The SIF inputs ir.'>y be stripped out of that SIF tape 
and sent over the data link to the Central Computer Complex. The operator 
must direct RTAT to submit a job to the data link. This is accomplished as 
f allows: 

(a) On the System Control Pcuiel, select the ON position for DAU 
DATA, set the Data Ident thumbwheels to 50, and press the ENTER button. 

RTAT will write end of file marks on the SIF tape, rewind it, read the input 
data from it, and leave the SIF tape positioned to continue writing on it. 
RTAT will send over the data link a batch job consisting of a standard 
set of control cards with the input data as a SIF. This set of control 
cards will convert the Sigma 3 core image SIF to a Control Data core image 
SIF and save it on the Central Computer Complex with a unique file name. 
Normal Central Conputer Complex jobs may then be sutxnitted to access this 
data file and process it. 
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Link from Central Computer Complex 

The Data Ident thumbwheel setting of 60 has been reserved for use in 
conjunction with the implementation of the data link from the Central 
Computer Complex to the Sigma 3* 

Skip RTAT Idle Loop 

It may sometimes be desirable to direct RTAT to cease its idle loop 
calculations. Tnis will not only remove a heavy computational load on the 
Sigma 3 but may permit OAP setup changes to be made without worrying about 
their impact on the RTAT idle loop. This may be accomplished in the following 
way. 

(a) On the System Control Panel, select the ON position for Dau 
Data, set the Data Ident thumbwheels to 90, and press the Enter button. 

RTAT will skip all execution including the idle loop. RTAT will continue to 
check for new Data Ident values which would require execution to begin 
again. 


RTAT Error Messages 

RTAT issues error and warning messages when it detects a problem in 
executing the requested function. The error messages appear on the 
typewriter and the warning messages appear on the line printer. A list of 
these messages along with the indicated problem is presented below; 

” NEED NEW RTAT SETUP .”- This message appears if RTAT attenpts to execute 
without having received an input setup deck. 

” DATA ID XX INVALID .”- RTAT does not recognize data ident xx. It will 
not attempt to execute unless it recognizes the data identification code. 
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” AVGREC ERROR nnnnn ” or " DIGICO ERROR nrmnn .”" These messages appear 
if RTAT was returned error code nnnnn by the designated GAP support library 
routine. Refer to the documentation of the appropriate routine for the 
Precise error encountered. If either of these messages appear repeatedly, 
it is likely that the data base has been clobbered, although it is possible 
that an GAP or hardware malfunction has occurred. It is necessary to re- 
boot the RBM system and restart the GAP and the RTAT from scratch to recover 
from a clobbered data base. This should be done when the original problem 
which clobbered the data base has been corrected. 

” NG SV NAME FGR CH nn .”- This message appears if there is no entry in 
the GAP channel name table for scanivalve channeJ nn. RTAT uses the first 
two characters of the GAP name to generate the valve name portion of the 
port leading names. An unnamed scanivalve channel is ignored by RTAT. 

” NG INTERACTIGNS . This message appears when RTAT is asked to compute 
weight tare factors without a balance interaction deck in the setup constants. 

The warning messages which appear on the line printer are self- 
explcuiatory. 


System Shutdown 

When all of the data have been obtained, the shutdown of GAP and RTAT 
is accomplished as follows: 

(a) Perform an GAP analog system calibration 

(b) Perform an RTAT system calibrate computation 

(c) To close out the GAP raw data tape: on the System Control Panel, 

select the GN position for EGT, and press the Enter button. Repeat this 


step three to six times 
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(d) To close out the RTAT output tape: on the System Control Panel, 

set the Data Ident thumbwheels to 30 and press the Enter button, set the 
Data Ident thumbwheels to 31 and press the Enter button. Repeat this step. 

(e) Rewind the data tapes, remove them from the tape drives, remove 
the write permit rings, and verify that the tapes are labelled correctly. 

(f) Press the Copy button on the graphics terminal hard copy unit 
to obtain a copy of the last plot. 

(g) Remove the output from the line printer. 

(h) Power down the Sigma 3 and the graphics terminal. The Data 
Acquisition Unit power must be left on. 
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Figure H-1.- Data identification codes 
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RTAT MANUAL CHECKPOINT GUIDELINES 

In order to verify that RTAT is correctly computing the data, it is 
necessary to ma. ually go through the calculations for a checkpoint. To 
assist in verifying intermediate calculations, it is recommended that a 
data base checkpoint be recorded for a wind off zero and for a wind on data 
point. Using these two points, a manual calculation should be made which 
starts with all of the input and computes all of the output. The algorithms 
which RTAT uses have already been discussed in sufficient detail to permit this 
manual calculation. However, the manual calculation is by no means trivial. 

It is the purpose of this section to outline some of the steps in the 
calculation for the case of a model with a single balance at angle of attack. 

(a) For the wind off zero, use the engineering unit equations to 
calculate initial loads 

(b) For the wind on data, use the engineering unit equations to 
calculate uncorrected delta loads 

[fu] 

(c) For simplicity, assume second order balance interactions are 
negligible to calculate correct delta loads 

Cf] 


308 



APPENDIX I 


(d) Calculate correct total loads 


(e) Use the engineering unit equations to calculate indicated 




(£) Calculate sting bending angles 


'^s'°‘s''*'s 


(g) Calculate balemce attitude with respect to gravity 


4> ,CL ,ll> 

g g g 


(h) Calculate model attitude with respect to wind 




(1) Calculate weight tares 


[ftare] 
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(j) Calculate aerodynamic loads on balance 

[fbal] 

(k) Calculate uncorrected tunnel parameters 

(l) Correct tunnel parameters for blockage and jet boundary 

(m) Calculate base and chamber pressures and corrections 

(n) Calculate model axis coefficients 

[cma] 

(o) Calculate stability axis coefficients 

[csa] 

(p) In order to check that second order interactions have been 

correctly applied, it is convenient to assume that the components are correct 

and verify the uncorrected components. This eliminates the iterative 

calculation required if the uncorrected components are assumed and the 

components are verified. This procedure assumes that [p] and [p ] 

o 

as computed by RTAT are correct and then proceeds as follows: 
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[ft] = [p] + [f ] 
o 


[f2t] = [ft][ft]^”'^®P°®® 


[fUT] = [ci][ft] + [c2][f2t] 

[f2 ] = [f ][f ]Transpose 

o o o 

[fU ] = [ci][f ] + [c2][f2 ] 

o O O 

[fu] = [put] - [fu ] 

o 

The second order interactions have been correctly applied by RTAT if 
[fu] computed by this procedure agrees with the [fu] used by RTAT within 
the accuracy specified. 

(q) The calculation of the angles using the Euler rotations presents 
a tedious manual calculation. An efficient manual algorithm is presented 
below. This algorithm consists of an initialization phase, a computation 
loop which is repeated for each rotation in sequence, and a termination 
phase . 

Initialization 
RV(1) = 1 
RV(2) = 0 
RV(3) = 0 
RV(4) - 0 
RV{5) - 1 
RV(6) - 0 
RV{7) - 0 
RV(8) - 0 
RV(9) » 1 
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Computation Loop 

Let A = rotation angle in degrees 

If A is a yaw rotation, let 
N1 = 1 
N2 = 2 

If A is a pitch rotation, let 
N1 = 1 
N2 = 3 

If A is a roll rotation, let 
N1 = 2 
N2 = 3 

Now compute 

T = RV(Nl)*COS(A) - RV(N2)*SIN(A) 
RV(N2) = RV(N1)*SIN(A) + RV (N2) *COS (A) 
RV(Nl) = T 
N1 = N1 + 3 
N2 = N2 + 3 

T = RV(Nl)*COS(A) - RV(N2)*SIN(A) 
RV(N2) + RV(N1)*SIN(A) + RV (N2) *COS (A) 
RV(Nl) = T 
N1 = N1 + 3 
IJ2 = N2 + 3 

T = RV(Nl)*COS(A) - RV(N2)*SIN(a) 
RV(N2) = RV(N1)*SIN(A) + RV (N2) *COS (A) 
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Repeat this computation loop for each rotation angle in sequence. 


Termination 

The resulteuit roll, yaw, and pitch angles, exclusive of 
quadrant determination, are given by 





= - arcsin (-RV(2)) 


e - 
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